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The series expansion for the evolution of the correlation functions of a finite system of hard spheres 
is derived from direct integration of the solution of the Liouville equation, with minimal regularity 
assumptions on the density of the initial measure. The usual BBGKY hierarchy of equations is then 
recovered. A graphical language based on the notion of collision history originally introduced by 
Spohn is developed, as a useful tool for the description of the expansion and of the elimination of 
£NJ , degrees of freedom. 

o , 

CN \ 1. INTRODUCTION 



In his famous derivation of the Boltzmann equation [8], O. E. Lanford makes use of a series expansion for 
the time-evolved correlation functions of a classical finite system of hard spheres in a box. This expresses the 
n— points correlation function at time t as a sum of integral terms involving all the higher order correlation 
functions at time zero. The expansion is derived, though not rigorously, from iteration of the BBGKY 
hierarchy of integro-differential equations, and is considered as a "series solution" of its Cauchy problem. A 
rigorous validation of the hierarchy and of the series has been given years later by H. Spohn in an unpublished 
'■■^ | note [TH, and by R. Illner and M. Pulvirenti in Q (see also the book Q), using different methods. 

In both the previous papers an assumption on the initial measure is made to derive the BBGKY hierarchy, 
that is the continuity along trajectories of the hard spheres flow. However, there is no physical reason to 
expect such a regularity property to hold, and it is worthwhile to notice that the final series expansion makes 
perfectly sense without assuming it. In fact, Spohn observes at the end of his note, by a density argument, 
that the expansion can be extended to a more general class of measures having no continuity properties. 
On the other hand, the interpretation of the BBGKY hierarchy as a family of partial differential equations 
is not at all easy, nor standard in any case, since it relies on the nontrivial properties of the operator T t 
of the hard sphere dynamics. Hence, the series solution concept appears to be more appropriate for the 
description of the dynamics in terms of probability distributions, and one wonders whether it is possible to 
derive it without going through the usual hierarchy. The present paper is devoted to a derivation of the 
series expansion for the correlation functions, which is not based on the iteration of the BBGKY equations, 
and never requires continuity along trajectories. We rather construct a method of direct integration of the 
solution of the Liouville equation, that allows to establish the validity of the expansion in a sense even 
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, stronger than those obtained in the existing literature: the result holds for all times in a fixed full measure 

invariant subset of the phase space, exactly as it happens for the existence of the dynamics of the underlying 
system of particles. The hierarchy of integro-differential equations is then recovered by resummation of the 
series, without additional assumptions on the initial measure, thus strengthening an analogous result in Q. 

Let us recall the derivation of Lanford and state our main result in an informal way. Consider the vector 
of correlation functions p — {p n }n>i, where p n is defined over the phase space of n hard spheres of mass to 
and diameter a > in a box A. A point in this space is an ?i— tuple (x\, ■ ■ ■ ,x n ),Xj — (qj,pj), specifying 
position and momentum of the n particles. If N is the total number of particles, we set p n — for n > N. 
Then the BBGKY hierarchy for the evolution of p can be written 

^-p{t) = Hp{t) + Qp{t), (1.1) 



where 



(Hp) n (xi, - ■ ■ ,x n ,t) = |i? n , Pnj(xi, ■■ ■ ,x n ,t) (1.2) 

is the n— particles Liouville operator acting on p n (including the effects of elastic collisions) and the collision 
operator is denned by 

(Qp) n (xi, ■ ■ ■ ,X n ,t) — a 2 ^ J dpdww ■ Pn+l(%l, ■ • ■ ,x n , q 3 + aw, p. t) . (1.3) 
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Here p is integrated over all M 3 , and w runs over the unit sphere. 

If t — > Tt(x±, • • • , x n ) is the flow of the dynamics, define the translation along trajectories of a vector of 
functions / = {f n } n >i as 

(S(t)f) n (x u - ■ ■ ,x n ) = f n (T_ t (x lr ■ ■ ,x n )). (1.4) 
Then, integration and iteration of Equation (1.1) leads to the formal solution 

pit) = S(t)p(0) + jr[dt 1 [ 1 dt 2 ---[ ml dt m S(t - t^QSih - h) ■ ■ ■ QS(t m )p(0) . (1.5) 
m=1 Jo Jo Jo 

In this paper we analyze in detail the structure of Eq. il.5\) and prove that it holds, for all times in a full 
measure subset of the phase space, for any absolutely continuous measure with density symmetric in the 
particle labels, and bounded by an equilibrium-like distribution. The hierarchy can be obtained then, in 

a mild sense, by taking the derivative. No assumption of continuity is needed even for this last operation. 
We also allow the total number of particles N to be non fixed by the initial measure. The boundedness 
requirement is stronger than the necessary, and it is the same used by Lanford to control the convergence of 
the series in the Boltzmann-Grad limit. Here it is made to control easily through all the steps the integrals 
over momenta of the type (jl.3p . (I1.5[) . 

The main interest of the discussion is the method of the proof. For n — N Eq. (jl.5[> reduces to the 
evolution of the density function, that is the solution of the Liouville equation: 

p N (xi,--- ,x N ,t) = p N (T_ t (xi, ■ ■ ■ ,x N ),0) . (1.6) 

It is desirable that we can construct the series expansion for the p n from direct integration of (jl.6|) over all 
the phase space of N — n particles compatible with a fixed state (aq, ■ • • , x n ). We show that in fact this can 
be done by eliminating the degrees of freedom one by one. To achieve the integration of the single degree of 
freedom, it is important to understand the structure of the right hand side in (|1.5j) . This has been widely 
studied since the work of Lanford Q , see for instance Q or [l3[ ■ It results that the integrand function in 
the generic term of the formula, depends on the states assumed by certain clusters of particles following a 
fictitious evolution: this is constructed from the state (x\, ■ ■ ■ ,x n ) at time t, by suitably adding more and 
more particles as the time flows backwards. Following [13J, we shall call collision history such an evolution. 

The collision histories can be represented graphically in terms of special binary tree graphs. Therefore, 
a graphical picture of the series expansion (|1.5[) is obtained. This representation is our basic tool. In 
fact, it turns out that the integration of a degree of freedom itself can be translated in graphical language, 
through appropriate operations over tree graphs. The graphical rules corresponding to the elimination of a 
single degree of freedom, clarify how the various terms of the expansion for p n emerge from those for p n +i, 
thus considerably simplifying the presentation of the proof. The analytical operations corresponding to 
these rules, are nothing but a suitable partitioning of the integration domain, and convenient representation 
(change of variables) of the subsets of the partition. Nevertheless, in order to establish the graphical rules, it 
is also essential to prove that some classes of collision histories give a net null contribution to the integration 
of the degree of freedom: this is done again with the help of the tree graphs, by showing explicit one by one 
cancellations among the collision histories of these classes. 

The paper is organized as follows. In Section [5] we define the model, we introduce our notations and 
state our assumptions on the initial measure. In Section [3] we introduce the concept of collision history, as 
well as the graphical rules for its representation, and explain how to represent formula {TT5J in terms of the 
tree graphs. In Section 2] we present our main results, while in Section [5] we discuss the proof of the main 
theorem, establishing the above mentioned graphical integration rules, and applying them to the generic 
inductive step. In Section [6] we present the conclusions and make some comparison with existing literature. 
Some technical aspects of the proof are deferred to the Appendices. 

2. THE HARD SPHERE SYSTEM 

In this section we set model and notations, which we inherit essentially from (l2l |. and state some prelimi- 
nary result on the hard sphere dynamics (Section 12 Al) . In Section [2 Bl we introduce the class of measures we 
will work with. 
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A. Model and notations 

Let us consider a system of N hard spheres of equal mass m and of diameter a > in a box A C R 3 . 
Denote Xi = (qi,Pi) eAxl 3 the configuration of the i-th particle, i = 1, • • • , N. A is bounded and has a 
piecewise smooth elastically reflecting boundary <9A. 

Between collisions each particle moves on a straight line maintaining unchanged its velocity. In a collision 
of two hard spheres at positions <?j,<7j with uj — (qi — qj)/\qi — qj\ = (qi — qj)/a G S 2 and with incoming 
momenta p\,p'j (that means {p^ — p'j) • w < 0), the outgoing momenta Pi,Pj (with (pi — pj)-ui > 0) are given 

by 

Pi = Pi - w[w • (pi -Pj)] , 

P'j = Pj + ' (Pi ~ Pj)] > (2- 1 ) 

as a consequence of conservation of momentum and energy. Moreover, in a collision of a particle with 
momentum p' i with dA at a regular point q (there is only one point of contact between the wall and the 
sphere) the reflected outgoing momentum pi is given by 

p i =p' i -2h{q){n{q)-p' i ), (2.2) 

where n(q) is the inner unit vector normal at q to dA. It is easy to see that the collision transformations 
(|2.1j) and (I2.2j) preserve Lebesgue measure on K 3 x R 3 and R 3 respectively. 
We may introduce the n— particle phase space, n = 1, • • • , N, 

r„ = {(xx, ■•• ,x n ) G (A x R 3 )" | | ft - q| > a/2 for every 

q G <9A, |% -qj\ > a, i,j = !,-■■ ,n, i ^ j} . (2.3) 

A state of the system is given by a point in the whole phase space Tjv. 

Under few simple regularity assumptions on dA (see [2[ for the details) the dynamics determined by (|2.1j) . 
(|2.2p and the free flow has been shown to exist in More precisely, it has been shown that there 

exists a subset T* C T„ of full Lebesgue measure dx\ ■ ■ ■ dx n such that for all t G K and for every point 
(xi, • • • , x n ) G r* the flow of the n— particle dynamics 

t^T} n \ Xl ,--- ,x n ) er; (2.4) 

is well defined. For all £ the mapping (xi, • • • ,x n ) — > T^ n \xi, • ■ • ,x n ) is uniquely defined as an invertible 
transformation from T* to T* . Moreover, Lebesgue measure on T* is preserved by the flow, being preserved 
at each single collision transformation of the type (|2.1j) . (|2.2p The flow can be extended to be a measure 
preserving map over the whole T n : we refer to [14| for a detailed discussion on the measurability properties. 

The set T n \ T* , which is of null Lebesgue measure, can be defined as the subset of all the points of T n 
which evolved in time run into either (see for instance |2j, page 16): 

• a "multiple" collision, that is simultaneous contact of more than two hard spheres or simultaneous 
contact of two hard spheres with each other and at the same time with OA; 

• a grazing collision with the wall (h(q) ■ p[ = 0) or a grazing two-body collision {{p[ — p'j) • w = 0); 

• a collision of a particle with a "singular" point of dA; 

• infinitely many collisions in finite time. 

The flow through such situations is not determined. We shall refer to them as the "singular configurations" . 
Some examples in which a particle undergoes infinitely many collisions in a finite time are given in [2j. 

Let us mention here a fact related to the properties of the flow. Call T the subset of dT n collecting all 
the multiple collisions, the simultaneous collisions (more than one collision occurring at the same time), the 
grazing collisions (between particles or with the walls) and the collisions with non regular points of dA. 
Consider the collision surfaces 

®n,ij = {(^1) ' ' • ) x n) G dT n \ T s. t. \qt - qj\ 

and qj = q t + aw, (p 3 - Pi) ■ W > 0} , 
®n,ij = {(^1) ' ' ' 3 3>n) € dT n \ T S. t. \qi - qj\ 

and q 3 = qi + aw, (p 3 - pi) ■ w > 0} , 



= a 



= a 

(2.5) 
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and 

= u,-* +( r ) 

*n,i _) = {(si,-- - ,0 6 5r„\T s. t. | ft - g | = o/2 

for some regular g € 9A, and • n(q) > (<)0} . (2.6) 

We have a decomposition of the boundary 

ar„ = cf>+ u *" U *+ U *" U T . (2.7) 

The Lebesgue measure on r„ induces, through the flow, a measure da n on dT n (||, whose restrictions 
onto $^ - , $^ 4 are given respectively by 

da n i j = ±dxi ■ ■ ■ dxi • • ■ dxj-idxj+i ■ ■ ■ dx n dpjdwa 2 w ■ (pj — pi) , 

da^.i = ±dxi ■ ■ ■ dxi-idx l+ i ■ ■ ■ dx n dqdpiPi ■ h(q) , (2.8) 

where w is the unit vector pointing from qi to Qj, OA 9 q,qi = q + h(q), and dq is the measure over the 
surface <9A. The set T has null a measure, and the Lebesgue measure on L„ can be written as da n dt, t being 
the time of the last collision in dT n . In references Q and [H it is proved that the flow (|2.4[) is well defined on 
9r„, again almost everywhere, with respect to the measure da. This existence property of the hard sphere 
dynamics is important for the derivation of the BBGKY hierarchy of integro-differential equations, as we 
will discuss in Section |4 Al fsee also Q). 

We shall collect the above results, for easy recall in the future, in the following 

Proposition 1. (Existence of the dynamics) The setY n \T* L C T„ defined by the above list is a null Lebesgue 
measure subset. Moreover, its intersection with the boundary dT n is a null measure subset of dT n with respect 
to the measure da n . 

For the proof, we refer to 0, Theorem II. B. 2, page 19 (see also @ and Q). 

As in 12j, we do not identify ingoing and outgoing momenta, but we regard them as corresponding to 

distinct points in phase space, so that the flow is only piecewise continuous in t. Then, when necessary, 
we distinguish the limit from the future (+) and the limit from the past (— ) writing 

T}±'(xi,...,x n )= lim T}?'(xi,...,x n ) . (2.9) 

E->0 + 

We list some definitions that will be useful in what follows. 

T N - n (xi,- ■ ■ ,X n ) = {(Xn+1, ' • ' , Xn) S TjV-n | \qi — qj\ > a 

for i = 1, . . . , N and j = n + 1, • ■ • , N} , (2.10) 

for (xi, • • • , x n ) E r„; that is rjy_ n (a;i, • • • , x n ) is the set of the possible configurations of N — n particles 
when we have n other particles in [x\, ■ ■ ■ ,x n ). We call Qi(xi, ■ ■ ■ ,x„,p) the points w on the unit sphere 
surface such that the configuration [x±, ■ ■ ■ , x n , qi + aw,p) is compatible with the hard core exclusion and it 
does not run into a singular configuration at any time: 

Qi(xi, • • • ,x n ,p) = {w S S 2 | (xi, ■■ ■ ,x n ,q l + aw,p) G (2-11) 

for i = 1, . . . , n, (xi, ■ ■ ■ , x n ) G T n and p€l 3 . If 

Cli(xi, ■ ■ ■ ,x n ) = {w G S 2 | (xi, ■ ■ ■ ,x n ,qi + aw,p) G r n+1 Vp G M 3 } , (2-12) 

then fli(xi,--- ,x n ) \ fli(xi,--- ,x n ) is a set of null Lebesgue-induced measure on S 2 for almost all 
(x\, ■ ■ ■ ,x n ,p) G r„ x M 3 , by Proposition [T] Furthermore we define Qi+(xi, ■ ■ ■ ,x n ,p) (Oj_(a;i, • ■ • ,x n ,p)) 
the points of Qi(xi, ■ ■ ■ ,x n ,p) corresponding to outgoing (incoming) collisions: 

n i+ (a;i, • • • ,x n ,p) = {w G Cli(xi, ■ ■ ■ ,x„,p)) \ w ■ (p - pi) > 0} , 

Qi-(xx, ■ ■ ■ ,x n ,p) = {w G fii(xi, • • ■ ,x n ,p)) | w ■ {p-pi) < 0} , (2-13) 

and analogous definitions for fij_|_(xi, •• ■ , x n ), f2j_(a;i, • ■ • ,x n ). 
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The following subsets of T* will be used to describe the time evolution of correlation functions: 

rj^ 0) = rjy = /Cat = f N = T* N and, for n < TV : 
It (0) - {x n G r; I for any 1 < k < N - n, it is x n+k G r^ +fc 
for almost all (x n +i,- • • ,£«+fc) £ r fe (x„)} , 

r+ = {x n g r; | for any 1 < k < n - n and s g r, it is (T^fej, x n+u ■ ■ • , x n+k ) g r; +fc 

for almost all (x n+1 , ■ ■ ■ ,x n+k ) G r jt (ri" ) (x„))} 

= H rW(r«°)) , 

= k„ e r; s.t. (T s (n) (x„), qj (s) + aw,p) e K n+l for all j = 1, • • • ,n 
and almost all (s,p,w) £ I x I 3 x (T^ n) (x n )) .} , 
f„ = rtf|/C„. (2.14) 

In the definition of IC n we put qj(s) = {Ts n \x_ n )) qj . The first two definitions ensure also that 
(T^ n \x n ),x n+ i 7 ■ ■ ■ ,x n+k ) G r| i+fe for every fc, s and almost all (x n+ i, • • • ,x n+k ) G T k (T^ n \x n )). Though 
it is not clear whether the two sets and T* coincide for n < N, we shall prove, as an extension of the 
result in Proposition [T] on the existence of the dynamics, that 

|r„\rt|=o, (2.15) 

where | • | denotes Lebesgue measure: see Appendix [5] where it is proved also that the restriction of the 
same set to dT n is da— null, and that it is 

|r n \/C„|=0. (2.16) 

From now on time t is always supposed to be positive, without loss of generality. We will use the short 
notation x n = x\, ■ ■ ■ ,x n and, when there is no risk of confusion, and we will simply call "particle z" a 
particle whose configuration is labelled by an index i. We shall set m = 1, since the role of the mass is trivial 
in all the discussion - see formulas (11.31) , (11.51) . 



B. Measures over the phase space 

Since all the particles of the system are identical, we will work with the space Cn of measurable functions 
In ■ Tat — ► R, symmetric in the particle labels (/jv(II(xi, . . . , x^)) — /at(xi, . . . , xn) for any permutation 
II), and having a boundedncss property of the type 

N 

\Mx 1 ,...,x ff )\<A'[[hp(p j ) , (2.17) 

3=1 

on Tat, for some A, ft > 0. Suppose to have an initial measure P on Tjv with density /at G Cn with respect 
to Lebesgue measure dx\ . . . dx at, 

P(dx\ . . . dxisr) — Jn(x\, . . . , XN)dx\ . . . dxN ■ (2.18) 

Then, because the flow T^ N ^ preserves the Lebesgue measure, the evolved measure at time t has a density 
f N (t) given by 

f N (xi,...,x N ,t) = f N (T^l(x u . . . , x N )) (2.19) 

almost everywhere in Fat, which is the Liouville equation in a mild form (the + sign is a convention). Points 
of Pat \ are removed from (|2.19[) . Estimate (|2.17p is preserved by the flow by conservation of energy. 
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Hence, /jv(i) G £jv- Of course since the flow is only defined almost surely, even densities that are 

regular at time zero will only be Cj^— functions at time t. 
We define the correlation functions p n , n = 1, 2, . . . by 



Pn{xi, ■ ■ ■ ,%n,t) = N. . . (N - n .+ f) / dx n+1 . . .dxNfjsrixx, . . . ,x N ,t) , n<N, 

JT N -n(xi,...,Xn) 

Pn = 0, n>N , (2.20) 

Pn j * * * j ^n) — Pn j * • • ) •> 0) 7 

where equality is in the space £ n , and points of T n \ Tn are excluded. Observe that 

n 

\p n (x 1 ,...,x n ,t)\<A , '[[h p {pj) , (2.21) 
3=1 

where A' can be taken equal to a pure constant times N n \A\ N ~ n . The volume of the system |A| and the 
diameter of the spheres a will be kept fixed along the whole paper, and of course by the hard core exclusion 
N will be bounded by 3|A|/47ra 3 . 

Let us say once and for all that, as in this paper we work with densities of measures, all equalities will 
hold for any fixed version of /iv,/jv(t) and p n {t) in their equivalence class, and all statements will be true 
only almost everywhere in Tjv or in its subspaces r„. Of course, the subsets where the involved flows of the 
dynamics are not defined for any time must be always excluded. In particular we will show that, assuming 
(I2.19j) and (12.20)) to be valid over the full measure subsets r„, all derived formulas (and in particular the 
final expansion) are still valid in the same set. If we like things to be more definite we can always think to 
fix a version of Jn, fN(t) assigning, for instance, zero value on the null set T n \ T n . 

We remark that P can be, in general, any measure with density in C n . In the case P is a probability 
measure, the quantity 
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dxi ■ ■ ■ dx n p n (xi, ■ ■ ■ ,x n ,t) (2.22) 

w 



N ■ ■ ■ (N - n + 1) 

is the probability of finding particles 1, 2, • • • , n at time t in the Borel set W 6 r„ 



3. COLLISION HISTORIES 



In this section we analyze the structure of the expansion on the right hand side of (|1.5I) . This is given 
in general by a large variety of terms. In each of these terms the integrand function contains a time-zero 
correlation function evaluated in a configuration of particles which can be found by flowing backwards in 
time the configuration x n , and suitably adding new particles at the times t± , t% etcetera. The new particles 
appear in a collision configuration with one of the pre-existent particles. This describes a special (fictitious) 
evolution that will be called "collision history" , a name first used by Spohn in [l2[ ■ 

In order to have a clear picture of the many terms of the expansion, and of the configurations of particles 
involved in them, we shall establish rules for their graphical representation. In particular, we will show that 
Equation (jl.5[) can be written as a sum over a set of tree graphs. We will introduce the convenient class 
of trees in Section 13 Al a class of decorated trees to be associated to the collision histories (as explained in 
Section 13 B| , and a class of trees with less decorations corresponding to the terms of the expansion. We will 
give the rules for this correspondence in Section 13 CI where also an explicit formula will be given for the 
generic term of the expansion. 

We want to stress since the beginning that the collision history is not a real trajectory of the particle 
system, and the associated collisions are not a sequence of real collisions. The correspondence between 
collision histories and sequences of real collisions is only very indirect (|12||). 



A. A family of trees 



We begin by considering binary tree graphs with generic node as in Figure [TJ one segment crosses the 
node while the other segment is generated by the node. In all the diagrams time will always flow from right 
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to left along a horizontal axis. We shall agree to draw the trees in such a way that the root corresponds 
to time t while the endpoints correspond to time zero, and that one of the two above mentioned segments 
attached to the node is horizontal, while the other has a (meaningless) slope between and tt/2. We call 
line of a tree the straight segment which left extremum is its generating node (or the root of the tree) and 
which right extremum is one of the endpoints of the tree. 

No trees will be considered with two or more nodes corresponding to the same time. Call 

A(x n ;[0,t}) 

for 1 < n < N,x n £ £ R, the set of all such binary trees with a number of nodes m variable in 

(0, • ■ ■ , N — n), exactly m + 1 lines (and endpoints), and no decoration other than the following: 

1. a label x n attached to the root; 

2. for n > 1, a label j £ (1, ■ ■ • , n) attached to each node crossed by the line ending in the root of the 
tree. 

We avoid to add a label t to the root of the trees if no confusion arises. See Figure [5] for an example. We 
may set A(x n ; [0, t)) = for n > N. 

We can always think to order the nodes of the tree from left to right with an index k — 1,2, ... ,m, which 
we call ordering number of the node. We refer to the line generated in the k— th node as the k—th line, and 
we refer to the line ending in the root of the tree as the root line. When the k— th node is crossed by the 
root line we denote jk the label associated to it. 

Two trees will be considered equivalent if they can be superposed, together with their labels and without 
altering their topological structure neither the ordering of its nodes. Hence, even though the nodes of a tree 
are not associated to precise values of the time, they are ordered along the time axis. For given number of 
nodes m and forgetting about decorations, there will be m! different trees in A(x n ; [0, t}); each of these trees 
can be decorated with the j labels in n m ° different ways, mo being the number of nodes crossed by the root 
line. 



FIG. 1: Structure of the generic node of a tree: one of the two lines is generated in the node, representing a new 
particle appearing in the collision histories described by the tree. 




FIG. 2: Example of tree in A(x n ;[0,t]). It will appear in the expansion for the n— points correlation function 
(n < N — 5); here G (1, • ■ • , n). In the figure < ti < ti < t. 

Now take a tree V £ A(x n ; [0,t]). Sometimes we will also use the notation V — (x n ,S) to remember 
that we fixed the root configuration. We shall call collision history and indicate it by T> — (x n ,6), the tree 
obtained from V by adding a triple (tk,Pkt ^k) £ M x I 3 x S 2 to the k—th node for every k — 1, 2, . . . , m, 
where 

• tk is a time variable, so that t m +i = < t m < • • • < ti < to = t; 
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• pk is a momentum variable; 

• the unit vector variable Wk has some complicated constraint depending on the other labels attached 
to the tree, which is defined in Section f3B I 

Call also 

the space of all the collision histories obtained in this way from A(x„; [0, t}) for n < N, and set A(x n ; [0, t]) = 
for n> N. See Figure [3] for some example. 




FIG. 3: Collision histories in A(a; n ; [0, t]) with 0, 1 or 2 nodes. 

Observe that our definition of collision history corresponds to the original one given in [12j . In fact 
A(x n ; [0, t]) is in a one by one correspondence with the subset of 

|J (N x E x E 3 x S 2 ) m (3.1) 

0<m<JV-ri 

given by the collections (to, ji, . . . , j m , t\, . . . , t m ,pi, . . . ,p m , w\, . . . , w m ) with the above mentioned con- 
straints over times and unit vectors, and with the variables jk £ N defined by 

jk = label attached to the k— th node, if the k-th node is crossed by the root line; 

jk = n + q , if the k-th node is crossed by the q— th line , 1 < q < k — 1 . 

(3.2) 

Notice that 1 < jf. < n + k — 1. In the notations T> = (x n , S), V = (x n , 5), we can identify S and S with the 
corresponding collections of variables: 

<5 = (m,ji, . . . ,j m ,t\, . . . ,t m ,pi, . . . ,p m ,wi, . . . , w m ) , 

6= (m,ji,...,j m ) . (3.3) 

B. The fictitious evolution of particles 

Let us now construct an evolution of particles 

£v 

to associate to the history V = (x n , S). The root of a tree is labeled by our starting configuration representing 
n particles at time t. In general the root line will represent these n particles from time to time t, and the 
fc-th line the (n + A:)-th particle, k = 1, 2, • • • , from time to time t^. The k— th node represents a binary 
collision between the particles associated to the crossed line and the generated line ((n + k)— th particle), 
and the triple (tk,Pk, £&fe) specifies the time of collision and the momentum and position of particle (n + k) 
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colliding: the vector joining the two particles involved in the collision and pointing towards particle (n + k) 
will be awk- Finally, the extra label j% 6 (1, . . . , n) in the nodes crossed by the root line tells us with which 
particle occurs the collision with particle (n + k). 

Given i = 1, 2, • • • ,n + m (m = number of nodes), we call 

Xi(s;T>) 

the configuration of the i-th particle at time s in the evolution £p (and we will put Xi(V) = a?j(0; V) for short) 
and we define this configuration for < s < t by the following construction. Take the root configuration 
x n G r*, put Xi(t;T>) — x n , and evolve it backwards in time as if there were no other particles in the space 

up to time t\ if m > (that is with the flow T_X t , ), and up to time if m = 0. This defines x n (s;T>) 
for ti < s < t. At time t\ stop your (n— particle) system and add particle (n + 1) in a state x n +i(ti;T>) 
with momentum p\ and position at distance aw\ from particle jx, with w\ G fX^ (x n (ti; T>),pi): at fixed 
x n ,ti we will have either an incoming or an outgoing collision between particles ji and (n + 1), depending 
on the chosen values of pi,w\. Then evolve backwards in time particles (1, . . . , n + 1) as if there were no 
other particles in the space up to time t 2 < ii (with T^l^l 2+ ); notice that soon after t± particle ji in the 
evolution S v will deviate from its free motion if and only if pi , wi correspond to an outgoing collision. At 
time t2 stop the system and add particle (n + 2) as above with momentum p 2 and position at distance au>2 
from particle j2 (defined by (J3T2J ) , with W2 £ £lj 2 (x n+1 (t2]T>) , P2) ■ Later on evolve your (n + 2)— particle 
system backwards up to time £3 < t2, and so on up to the final step, which is the evolution of particles 
(1, . . . , n + m) with the flow t!™^™^ from time t m > to time 0. We stress that the configurations x n (s; V) 
are always constructed by taking limits from the future. An example is pictured in Figure |4j 




FIG. 4: Trajectory drawn by the particles in a collision history of the type {x n ,52,\) of Figure [3] in the case 
n = 2, ji = 2,w\ £ fl2-,j2 = 1,W2 6 fh+. 

In the following we will call 

£v{s) 

the configuration of all the particles of £t> at time s, without specifying the number of such particles, so 
that Ex> = {£v{s)}o<s<t- In particular if s coincides with the time t^ associated to a node, then £t>(s) is 
the configuration of the particles of the evolution after having added the new particle generated in the node: 
£v{tk) = (x n+k ^i(tk; / D),qj k (tk;'D) + awk,Pk)- We call 

N v (s) 

the number of particles of £p at time s, Nx>(s) G (n, . . . , n + rn), and we name cluster of particles of £x> the 
time-dependent collection of particles described by the evolution. 

If we would prefer formulas instead of trees we should write, referring for instance to the trees in Figure [3l 

x n (s; x n , 5 ) = T^ + (x n ) , t > s > , (3.4) 

x n (s;x n ,8i) = TrJ + (x n ) , t > s > -t + ti , 

x n+1 {x n ,8i) = Tr^l\Tr^ +tl+ {x n ),q h {ti]X n ,5-i) +aw 1 ,p 1 ) , 
x n+2 (x n ,82,i) 

= T -tt+ ( T - n tt+t2+ ( T -'i+i 1 +fen) 5 9ji(*i;Sn> 5 2,i) + awi,pi) ,qj 2 (tr,x n ,S 2 ,i) + aw 2 ,p 2 ) , 

and so on. 
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Remarks. 

(1) Clearly all this construction is not well denned for x n e T„\ T* ; definition (|2.11j) ensures that the 
added particles do not run the system in a singular configuration. 

(2) As we anticipated above, the evolution £x> is not a real trajectory although it is constructed with 
pieces of possible real trajectories, and a collision history T> is not a sequence of real collisions of the system. 
In particular, notice that the configuration £ n +fe (s; T>) given by T> is defined only for times < s < tk, and 
that in general 

£v{t k ) 

is different from the limits of £p(s) as s — > ±tk. 

(3) Given the graph (x n , 5), if particles jk and n + k are attached to a node with time index tk, then in Ed 
for times s <E [0,tk] they can collide many other times between them and with the other particles appearing 
in the graph at those times, but not with other particles of the system that do not appear in the graph. In 
general any two particles appearing in the graph at a given time can be in a collision configuration. 



C. A graphical representation of formula (11.51) 

Suppose V 6 A(x„; [0, x n € T*, to have m nodes, m 6 (0, ••• ,N — n), and order them from left 
to right with the index fc, 1 < k < m as in Section 13 Al The tree will be associated to collision histories 
V G A(x„; [0,i]), obtained attaching to the nodes of V triples (tk,Pk,u>k), as explained in the same section, 
and - see Eq. (|3 - 3[) the two graphs will correspond to collections of variables 

V = (x n ,m,ji,--- ,j m ) , 

V = (V,tl, ■ ■ ■ ,t m ,Pl, ■ ■ ■ ,Pm>l, • ■ • ,w m ) , (3.5) 

where ji , ■ • • , j m are defined via (|3.2[) . 
We denote 

V{V) 

the value of the tree T> given by the rules summarized in what follows. Going from left to right, i.e. climbing 
the tree: 

• associate to the k— th node of the tree a weight factor 

W k (V) = a 2 w k ■ { Pk -p jk (t k ;V)) ; (3.6) 

• associate to the fc— th node of the tree an integration over a subset of M x IR 3 x S 2 given by 

dw k (3.7) 

, +fe _i(*k;i > ),Pfc) 

(remember that t = t$) where diik is the natural induced measure on f2j fc ; 
and at the end 

• associate to the m + 1 endpoints of the tree the (n + m)— th correlation function at time zero evaluated 
in the final configuration of £x>, that is 

pn+ m {xi{V), ■■■ , x n+m (V)) . (3.8) 

Hence V(T>) can be seen as an integral over times, momenta and unit vector variables attached to the 
nodes of the collision history T> £ A(x n ; [0, t]), of a product of a correlation function times a weight function 

m 

W(V) = ]J W k (D) . (3.9) 
fc=i 









dpk \ 




r 3 Jn 
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Explicitly, 

V(V) = f dh [ dp! f <Mi f dt 2 f dp 2 f dw 2 ... 

JO Jr 3 J^j 1 (x n (t 1 ;V),p 1 ) JO Jr 3 J^i 2 {x n+1 (t2\V),p 2 ) 

dt m / dp m / dw m R(D) , (3.10) 

o Jr 3 Jn jm {x 7l+m _ 1 (t m ;V),p m ) 

where R is called value of a collision history and is denned by 

R{V) = W(V) 

Pn-\-7n (^1 (©),.•• ,x n+m (V)) (3.11) 

which will be, in our assumptions (see Section [5]), a measurable function over the domain of integration in 
(13.1 Op for almost all x n £ T n . Formula p. 101) is a representation for the generic term of the expansion (jl.5|) . 
The domain of integration in f|3 . 10[) is maximal for all times if x n G K, n , see f|2 . 14|) . In fact K, n is defined 
as the maximal set of values of x n for which the evolutions Z^ x ^ (hence their value R(x n ,S)) appearing 
in the evaluation of V(V) are well defined for almost all values of S, that is for almost all times, momenta 
and unit vectors associated to the nodes of the evolution, and compatible with the hard core exclusion. The 
corresponding integrals in d5 are then extended over full measure regions over the sets compatible with the 
condition of hard core exclusion, so that we can substitute (• • • ) with Q, j i (■ ■ ■) in the expression (|3.10l) . 

The so defined V(T>) and R(V) can be seen as operators respectively over the spaces of variables 
{A(x n ;[0,t}) | n = 1,2,--- ,x n E T*,t > 0} and {A(x n ;[0,t}) \ n = 1,2,--- ,x n £ T* n ,t > 0}, wi th va lues 
in some space of functions over T* x (0, oo), n = 1, 2, ■ • • Notice that the integrals in £j and Wi in (|3. 10[) are 
over finite regions, while, in the assumptions of Section^ the integrals over pi are controlled by the estimate 
\p n + m (xi(V), ■ ■ ■ ,x n+m (T>))\ < (const) YTj=i hpiPj) Iljli hp(Pj)i which follows from (|2.21[) and conserva- 
tion of energy, and |W^(D)| < (2a 2 ^/^™ =1 p? + T^ = iP]) m ■ Hence the integrals in (13. 10)) are absolutely 

convergent, they define a measurable function over T n for any t > and, using \p\ q hp(p) / hp/ (p) < const for 
q > 0, &' < P, we have 

n 

\V{V)\ <ClIV(Pj) - (3-12) 

i=i 

with j3' < (3 and C depending on N, a, A, t. Finally, oserve that V(T>) is symmetric for exchange of particles 
Xi <H» Xj,i,j G (].,••• ,n), and simultaneous change i — > j,j — >• % in the value of the node labels, so that 
certainly the sum over all trees X^dgAO -[o t\) V^P) ^ s m 

In our main theorem we will prove that, starting with an initial density /at G Cn, this sum does give the 
time evolution of the correlation functions. 



4. THE EVOLUTION OF CORRELATION FUNCTIONS 



In what follows we present our main theorem. After the statement of the theorem and some general 
comment, we derive the usual BBGKY hierarchy of equations (Section 14 A[) . Finally, we present also an 
extension of the result to measures of grand canonical type (Section 14 Bp . 

Theorem 1. Let P be an initial measure with density /jv G Cn ■ Then for any t > 0, the time-evolved 
correlation functions are given by 

Pn(x n ,t)= V <P^> UeN > 

V^K{x n -[0,i\) 

almost everywhere in T n . For any chosen version of f^it), p n (t) satisfying \2. 1 9\) and \2. 20\) over the whole 
sets r„, the expansion holds for all x n £ T n and t > 0. 

In reference [12| this formula is called the "time-integrated form of the BBGKY hierarchy". Actually it is 
the complete expansion of the n— th correlation function at time t in terms of the higher order (n + m, m > 0) 
correlation functions at time zero. The number of terms in the sum is of course finite. 



12 



Remembering what has been said next to (|3.1|) . we may define a measure dS on A(x n ; [0, t)) as the counting 
measure with respect to the discrete variables m,jx,... , j m , and the Lebesgue measure with respect to the 
variables t±, . . . , t m ,pi, . . . ,p m , u>\, . . . , w m . With these notations we can say that R{V) — R(x n ,S) is a 
d5— summable function on A.(x n ; [0,t]) for all t and almost all x n G T n , and we can rewrite (|4.1I) as an 
integral over collision histories: 



almost surely in Y n . 
Remarks. 

(1) We do not need /jv and p„ to be continuous along trajectories of T t , that is we do not need 



for a. a. x N G Tat, where both the limits from the future and the past are understood. In particular, our 
initial density can distinguish between pre-collisional and post-collisional configurations (which can not be 
true if we assume for instance (|4.3[) on all T*). It is easy to show (see [12J) that, if the continuity along 
trajectories is assumed to be valid for /at, then the Liouvillc Equation (|2.19p together with some integrable 
bound on /jv imply that: (i) the same continuity property is also valid for /jv(i) and for p n (t) at any time 
t > 0; (ii) for almost all x n the map t — > p n (x n , t) is continuous. All these properties, even if assumed, would 
be not helpful in the discussions of the present paper. In [ljj and in [f| the continuity along trajectories is 
used to derive the series expansion (|4.1I) . 

(2) The bound (13.12[) follows from rough estimate of the right hand side of p.lOj) . as already explained. 
That is sufficient for our purposes. From the proof of the theorem it will be clear that the bound (|2.17l) could 
even be substituted with a weaker one, since it is just needed to ensure absolute convergence of the integrals 
in (|3.10j) . Our choice of the decay behavior for high momenta is the same used in the careful estimate of Q 
of the right hand side of (|3.10j) (see the details in Q), necessary to perform the Boltzmann-Grad limit: if 
the correlation functions satisfy \p n (x n )\ < c(Nz) n YYj=i hpiPj) for some c,z,/3 > 0, then the right hand side 



of g3D is bounded by | £ ro > Eu e x<»0 (B ;[ o, t] ) V ^ ^ c '( Nz T n" =1 V(Pi) £ m >o( const ■ iVa 2 ^)™ for 



some c',/3' < /3 and z' > z(/3'//3) 3 / 2 (here A m (x n ; [0,t]) is the subset of trees with m nodes). This ensures 
convergence for N — > oo, No? fixed, at least for sufficiently small t. 

(3) Our result is actually stronger than the one obtained in [l2[ via density arguments, which is the same 
expansion integrated over every Borel set in T n (and corresponds to the first statement in our theorem) . We 
know that there exists a full measure subset of the phase space where the dynamics of the hard sphere system 
exists for all times (Proposition!!]). Theorem[T]recovers this property for the evolution of correlation functions: 
the expansion (|4.1I) is valid for all times in T n , that is a full measure subset of T„, and invariant under the 
flow. This subset - see the definition of r„ in Eq. (12. 14)) - has not been characterized in a constructive 
manner: this would depend on details of the dynamics that have not been investigated. However, it will be 
clear from the proof that l\ is the maximal subset of the phase space where the result can be derived for all 
times. In particular, the second statement of our theorem is still true if we replace T n with any full measure 
invariant subset of it, say T-L ni satisfying the following "chain property": if x n 6 H n , then (x n ,y ,) G Hn+k 
for almost all y G Fk(x n ). 

(4) Choose a version of /at(^), p n {t) satisfying (|2.19|) and (|2.20p in a set H n C T n (it is sufficient H n C r£) 
as described in the previous remark. Call (remember definitions (|2.5|) . (|2.6I) and (12.71) ') 




n G N 



(4.2) 



lim /Ar(T s w (xi, • • • ,a;jv)) = /jv(^i, ■ • • ,xn) 



(4.3) 




(4.4) 




(4.5) 



for all times t > 0. The converse is also true. Here the restriction to Ti. n corresponds to the conventional 
+ sign used in writing the Liouville Equation (|2.19l) . In the next section we will use the above formula as a 
starting point to derive the integro-differential BBGKY equations. 
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(5) Theorem Q] and Lemma [2] of Appendix lAl immediately imply that, for any chosen version satisfying 

(|4.5p in Tn, the identity (|4.1I) holds for every time over almost all dT n . 

(6) Formulas (j4.ip and (|4.2j) suggest an interpretation of the contribution of a collision history to the right 
hand side in terms of constructive or destructive correlation effects of the "external" particles (i.e. those 
different from (1, • • • , n)) on particles (1, • • • , n) during the time interval [0, t]. Consider for instance the 
history (x n ,5{) = (x n ,l,ji,ti,pi,u)i) of Figure [3j This gives a positive or negative contribution to the right 
hand side of (|4.1|) depending on the sign of the weight factor Wi(x n) 8) associated to its node. In the first 
case, Wi(x n ,8) > 0, the (only) external particle n + 1 appears in an outgoing collision with particle j\\ its 
effect on particles (1, • • • , n) is that of creating the configuration (n, . . . ,x n ) at time t, in the sense that 
if we forgot the interaction effect of particle n + 1 on particle j\ at time t\ (thus modifying the trajectory 
drawn by the collision history), then by evolving forward in time we would not get . . . , x n ) at time t. In 
the other case, Wi(x n ,S) < 0, the particle n + 1 appears in an ingoing collision with particle j±: its effect 
on particles (1, • • • , n) is that of annihilating the configuration (x\, . . . ,x n ) at time t, in the sense that if 
we took into account the interaction effect of particle n+1 on particle j\ at time t\ (thus modifying the 
trajectory drawn by the collision history), then by evolving forward in time we would not get {x\, . . . ,x n ) at 
time t. More generally, in any tree, we can say that a node with positive weight factor describes a collision 
that creates, in the ordinary verse of time (that is going towards the root) a particle entering the next node 
of the tree - or entering the root if the node is the last one (i.e. creating x n at time t); while a node with 
negative weight factor describes the annihilation of such a particle. It is then clear why, for instance, in trees 
with two nodes, two annihilation weight factors (negative) correspond to a net positive contribution to the 
right hand side of (14.1[) . two weight factors of different type (one positive and the other negative) to a net 
negative contribution, and so on. 



A. The BBGKY hierarchy 

We want to show here how the usual BBGKY hierarchy of integro-differential equations is recovered from 
the expansion (|4.ip . We present below the bulk of this derivation and refer to the appendices for some 
technical details. We begin by fixing a version of the density and the correlation functions satisfying (|2.19[) 
and (|2.20p . for simplicity, on the whole f„, so that Eq. (|4.5[) holds over all fn and for all i > (everything 
that follows would hold also replacing T n with any subset T-L n : see Remarks 3 and 4 in the previous section) . 
As we will see, starting from that formula it is easy to obtain informations about the function of time 

* — ► Pn(Tf U \xi, ■ ■ ■ ,x n ),t) , (4.6) 

and its derivative, without additional assumptions on the initial measure. To begin with, it can be shown 
that, for any n and all x n £ fn~^ (or r^ + ^), the function is continuous for every t > - see Appendix |D| 

Secondly, we can rewrite the expansion (|4.ip in a resummed form, which is convenient to obtain infor- 
mations about the derivative, as explained in the following. Fix n < N, and rewrite the expansion (|4.1I) 
as 

Pn{x n ,t) = p n (T^ + (xi, ■ ■ ■ ,x n )) + / dti / dpi / dwia 2 wx ■ (px - p h (tx;x n ,8-i)) 

ji = 1 J Jrs JO jl (TW (i + ( S „),y 1 ) 

■Pn+i (T% i+ (x n ),q jl (t 1 ;x n ,S 1 ) + at&x.px)) + £ V(V) , (4.7) 

veA(x n ,[o,t]) 

m(D)>l 

for x n 6 r„, where (x nl 8\) is the tree with one node in Figure [31 and m(T>) is the number of nodes of T>. 
Remind that we can always substitute Q,j t with in the above expression (see Remark 3 in the previous 
section). In formula (I4.7P we wrote explicitly the lowest order terms (zero-nodes and one-node trees) of the 
expansion. Since we restrict to x n G K, n , in the integrals corresponding to the one-node trees we may use 
again Equation (|4.ip to substitute p n _|_i(-,0) with p n+ i(-,ti) : it follows that the second term in the right 
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hand side of (|4.7[) is equal to 



pt p p 

V] / dt x \ dpi I dw^wx ■ (p l -p^ihix^Si)) (4.8) 



ji = 1 J0 JR3 JSl^r.^^JA) 

■Pn+i (T ( ^ +ti+ {x n ),q n (tux n ,5i) + awx.px.tx 



V" / dt i d Pi / dwia 2 wi ■ (pi -p^fax^Si)) V] 

Ji = 1 -t+ti + V-^/'^w UeA(x ,,(ti;x 



V(V) 

CeS^ n+1 (ti;x n ,<5i)j[0,ti]) 
m(Z?)>0 



The second line of the last formula gives all the trees with at least two nodes, i.e. — ^(^)- Hence 

veA(x n -[o,t]) 

m(V)>l 

for x n G f„ we have found 

" pt p p 

Pn{x n ,t) — p n (ri"+(xi, ■ ■ ■ ,X n )) + / <f*l / #1 / ^1 

il= l-/0 JR 3 ^(T^^xJ.pi) 

•a 2 wi • (pi -p J - 1 (ti;^„,5i))p„+i ^ri"^ tl+ (x„),g J ' 1 (ti;x„,5i) + auJi,pi,ti^ (4.9) 

(which is again an absolutely convergent integral). Formula (|4.9|) is the resummed form of the expansion for 
the correlation functions, in the sense that iterating the equation N — n times we are back to the Equation 

dm. 

Recalling the continuity property stated at the beginning of the section, we can write also 

^ U )lnn \ -A — n Inn \ I / M. (Hi . . o . . \ lT^ n ) 



Pn(Tr ($J,t) = Pn(x n ) + / dh (Q n+ ip n+ l) (T™(jBj,h) (4.10) 

JO 

for x n G Tn i where the collision operator Q n +i acting on the time-evolved correlation function, is defined 

by 

(Qn+iPn+i) (x n ,t) = y^a 2 / dpi dww ■ (p - pj) p n+1 (x n ,qj + aw,p,t) . (4.11) 
j= i "'R 3 Jn i (x 1 J 

over r n x [0, oo). The integrand in (|4.10l) is a measurable function in the variable t\ for all x n G /C„, while, 
for all t, the definition (|4. 1 1[) can be extended to r n , providing a function in the space C n . The definition 
is, in our assumptions, independent on the chosen version in the sense that, if p n +i{x n+1 ) = Pn+i(x n+1 ) 
for almost all x n+1 G r n +i 5 then the same is true for the time-evolved functions for all £ > (see (|4.5j> ) 
and, by the continuity property stated in the second part of Lemma [5l the two functions coincide also 
for almost all (x n+1 ,t) G dT n+ i x [0, oo) (see also the final paragraph in the proof of Lemma [2]), so that 
Qn+ipn+i = Qn+iPn+i for almost all (x n ,t) G T n x [0,oo). 

Formula (|4.10l) shows that t — > p n (T^ (x n ),t) is also absolutely continuous. As a conclusion, we can 
state 

Corollary 1. Given an initial measure with density G Cn, and /jv(t)> Pn(t) satisfying &2.19\) and \2. 2(A) 
on f n , the function t — > {Q n +iPn+i) (T^ (x n ),t) is measurable and the correlation functions satisfy 

j t Pn{T^\x n %t) = {Q n+lPn+1 ){T^\xJ,t), neN, (4.12) 

for all x n G f „ and almost all t > 0. 



The subsets of the phase space involved in the assertion of the lemma have full Lebesgue measure. Remind 
lat fi +) is 
Remarks. 



that r^ +) is mapped by T t (n) onto T n for t > 0. 
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(1) If, additionally, Jn £ C(T n ) then, for all x n G f „ \ dT n and almost all t > 0, it could be proven that 
the right hand side is continuous in t. The boundary dT n is discarded as it contains the (possible) points of 
discontinuity along trajectories of the time-zero correlation functions. 

(2) We did not used the continuity along trajectories of the correlation functions. Thus the result strengthens 
the analogous in Q . Weaker versions of the hierarchy have been already proved without the assumption of 
continuity along trajectories, see @ or Q- 

(3) Unlike the series solution (I4.1[) . the BBGKY hierarchy in differential form explicitly involves restrictions 
of the correlation functions to sets of codimension 1: this has made crucial the property of existence of the 
flow on the collision surfaces (second statement in Proposition [1]). 



B. Indefinite number of particles 

Finally, we can extend the result of Theorem [I] to a more general class of measures with non definite (but 
finite) number of particles. We follow [l2[ for this purpose. Consider the grand canonical phase space 

r = u„> r„. (4.13) 

Then it will be r„ = for n larger then [3|A|/47ra 3 ], because of the hard core exclusion. 

Call C the space of measurable functions / : T — »• R, / = {fn}^Lo, symmetric in the particle labels 
(/ n (n(xi, . . . , x n )) = f n (x\, . . . , x n ) VYi, for any permutation IT), and having the boundedness property 

n 

!/„(£!,..., Ol < Ajlizhpfa)) , (4.14) 
i=i 

on r„, for some A,z,(3 > 0. We can put /„ = for n > [3|A|/47ra 3 ]. If P denotes a measure on T with 
density / G £ with respect to Lebesgue measure, then the time-evolved measure at time t has a density 
f(t) G C given by 

f n { Xl ,...,x n ,t) = f n {T^ + {x u ...,x n )) , neN (4.15) 

almost everywhere in T n . 

Given / E £, we define the correlation function vector p : V — > M, p — {pn}^Lo by 

°° 1 f 

p„(xx,...,x n ,t)=y^ — / dx n+1 ...dx n+k f n+k (x 1 ,...,x n+k ,t) , (4-16) 



where equality is in the space C n . Again we have that p e C and furthermore, the map defined by (14.161) 
has the inverse 



°° (— l) fc f 

f n (x\,...,x n ,t) = / dx n+1 ...dx n+k p n+k (x 1 ,...,x n+k ,t) . (4.17) 

k=Q JT h {xi,...,x n ) 



The following extension will be an immediate consequence of the analysis developed in the next section. 

Corollary 2. Let P be an initial measure on V with density f 6 C. Then for any t > 0, the time-evolved 
correlation functions are given by 

Pn(x n ,t)= J2 V( P), neN, 
veS(x n ,[o,t]) 

almost everywhere in T n . For any chosen version of f n (t), p n {t) satisfying J^._?5[ ) and {4-16ty over the whole 
sets r„, the expansion holds for all x n G T n and t > 0. 

(Here f„ is defined as in (|2.14[) with k > 1.) 

Each term in the sum in Equation ()4.16|) may be dealed with the procedure explained in Section [5] This 
leads directly to a tree expansion of the type in the right hand side of (14. 1| . in which the value of the tree, 
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say V(V), is computed in a slightly different way in terms of the density function /. To evaluate V(T>), 
follow the rules introduced in Section f3 CI substituting formula (|3.8j) with 

■ ■ ■ dXn+k 

(it - my. Jr k _ m ( Xl (T>),- ,x n+m (v)) 

■fn+k (xi(D), ■ ■ ■ ,X n+m (p),X n +m+l, • ■ • ,X n +k) • ( 4 -18) 

Performing the sum over k, that is Y^k>mi an< ^ usm S (HUH), we obtain the corollary. 



5. PROOF OF THEOREM |TJ 



In this section we prove our main result. We shall proceed by induction on n: supposing the statement of 
the theorem true for the function p n +i, we derive the expansion for the p n by integrating a single degree of 
freedom. The rigorous integration procedure is rather technical but, in spite of lengthiness of formulas (for 
which sometimes we refer to the appendices), the integration of a degree of freedom in a single term (tree) 
of the expansion admits a quite simple graphical representation in terms of "extraction" of subtrees and 
"reattachment" of extracted subtrees. These operations over trees are introduced in Section |5 A[ while the 
graphical integration rules are summarized in Proposition[5]in Section \5 Bl Along the proof of the proposition, 
in the same section, we present the analytical operations depicted by the operations over trees: they consist 
essentially in appropriate partitioning of the integration domain and representation of its subsets. However, 
this is not sufficient: to prove the proposition it is also essential to notice that a certain class of collision 
histories gives a net null contribution to the integral, because of one by one cancellations. This will be done 
in Lemma [T] Finally, in Section 15 CI we conclude the proof of the theorem, by discussing the summation of 
all the graphical terms obtained through the integration procedure. 



A. Tools: manipulation of trees 



The integration of degrees of freedom will be described by a manipulation of the trees involving "pruning" , 
"extraction" and "growth" operations, for which we will need some more notations. As in the previous 
sections, we will indicate with m = m{T>) = m(5) (or m(T>)) the number of nodes of the tree T> = (x n ,S) 
(collision history T>), and with jk — jk(T>) = jk{5) (jh(T>)) the variable defined in (j3.2p . Moreover, we will call 
A n , m the set of trees in A(x n ; [0,t]) with m nodes and not specified time and initial configuration (no labels 
attached to the root); clearly an element Q £ A„ jm is identified with a set of variables (n, m, jx, ■ ■ ■ ,j' m ), see 
(|3.3l) . For the trivial tree, that is the only one belonging to A^q, we will use the symbol T. 

Given V £ A(x n ; [0,i]), order its nodes from left to right with the index k as in Section [3l We name 
T>k £ Ai ,„' the subtree generated in the node number k, and we call T>/u £ A(x n ; [0, i]) the tree obtained 
from T> by pruning T>k (it will be m' < m{T>) — 1). 

Now we extend these definitions to the case k = 0. We call £>o:j,j = lj • ■ • ,n the tree in Ai iT7l / ob- 
tained from T> £ A(x„; [0,t]) by pruning all the subtrees T>k such that the node k lies in the root line 
and has a label jk j (it will be m! < m{T))). Similarly, we call D/ .j,j = 1, ••■ ,n the tree in 
A(xi,-- - , xj-i, xj+i, ■ ■ ■ ,x n ;[Q,t]) obtained from V £ A(x„;[0,i]) by pruning all the subtrees Vk such 
that the node k lies in the root line and has a label jk — j- Notice that in the case there is no node label 
with value j, it is ml = 0, Z>Q;j = T, and "D/o-j is the tree in A(xi,--- , Xj-\, Xj+i, • • • ,x„;[0,t]) which 
is identical to V except for the label attached to the root, i.e. 2?/o ; j — (^lr - ' > x j—ii x j+ii'" i x n ,5) if 
V=(x n ,5). 

We can also visualize the trees X>Q;j in the following manner. Imagine that the root line oil) £ A(x n ; [0, t]) 
is composed by n coincident identical lines, numbered from 1 to n and associated to the particles of the 
initial configuration x n . The j—th of these line, j — 1, • • • ,n, is thought as attached only to the subtrees 
generated in the nodes of the root line carrying a node label with value j. Then we can say that the subtree 
T>o.j, is obtained by extraction of the j—th line, together with the subtrees attached to it, from the tree T> 
(and by deleting decorations). The j—th line will become the root line of the extracted tree. What is left 
of the original T> after the extraction of T>o,j and after deleting the root label Xj as well as the node labels 
with value j, is exactly the tree "D/ .j, 



17 




FIG. 5: Notations for subtrees, pruned trees and extracted subtrees. 
Finally, for n < N, we define the composition of trees o^, k — (fco, • • • , k q ) <E Z 9+1 , 1 < fco < fci < ■ ■ ■ < k qi 

by 

%,-,*,;i:%„;[0,i])xA 1 ,,4A(x„;[0,t]) 
T> °k,i G = H such that 1-L ka = G, H/ ko = T> , 

the nodes of V.k have ordering 

numbers fei, • ■ • , k q in %, and 

the label jk of Ti is equal to i , 

for k q < to (25) + g + 1 and 1 < i < n + k - 1 , (5.1) 

and otherwise. This means simply that grows the tree V by attaching to it the tree G in such a way 
that: 

1. the root of G is attached to the root line of V when i E (l,-- - ,n), and to the r~ th line when 
i = n + r, 1 < r < m(D); 

2. a node with ordering number fc is created in the previous operation; 

3. the ordering numbers of the nodes of G in the resulting tree are given by (from left to right) 
fci, • • • ,k q ,q = m(G). 




FIG. 6: Examples of composition of trees: in V o^ ;i Q the indices fco, i indicate to which line of V (and between which 
nodes) has to be attached the root of Q, and the fci, • • • , k q indicate how to order the q nodes of the subtree Q in the 
resulting tree. 
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B. Integrating a degree of freedom 
Formula (|4.1|) is trivial for n > N, while for n = N it gives, graphically, 



PN(x N ,t) = x N = Pn(T^I(x n )) (5.2) 

a.e. in Tat, which is implied by (|2.20|) and (j2.19[) . 

We shall proceed by induction on n to prove the theorem: from (|2.20[) follows 

Pn{x n ,t) = -rr — / dx n+1 p n+1 (x n ,x n+1 ,t) , 1 < n < N , (5.3) 

N - n Jr 1 (x n ) 

so that assuming ()4.1|) valid for p n +i we can write 

Pn {x n ,t) = ^— [ dx n+1 V( V) (5-4) 

a.e. in r„, or exactly in F„ if we are assuming (|2.19l) and (|2.20p to hold over it: this will be understood in 
what follows from now on. Then we need to explain what is the result when we integrate the single degree 
of freedom in a tree of the set A(x n+1 ; [0,i]), i.e. we have to compute 



I(V) = I(V)(x n ,t) := dx n+1 V(x n+1 ,S), D = (x n+1 ,8)eA(x n+1 ;[0,t}) (5.5) 

for a set of x n of full measure in T n . Notice that the integral in the above formula is well defined as a 
measurable function over r„, see the final comments in Section [3 CI 

The computation of (I5.5[) will be the main part of the proof, and the rest of this section. After that, we 
must just sum the result over all the trees of the family A(x n+1 ; [0,/;]). 



1. Integration of a degree of freedom in a single tree 
Given V G A(x n+1 ; [0, f]), and selected a particle j € (1, • ■ • , n + 1 + m(V)), call 

q[ J \q { 2 3 \--- (5.6) 

the ordering number of the nodes, in T>, that belong also to the subtree with root line given by the line 
associated to particle j: 2?o;j in the case j € (1, • • • , n + 1), or 2\- in the case j = n + 1 + k, k > 0. 
Write gC?) = (q[ 3 \q^\- ■ ■ ). Define a variable I* = l*(V) by 

I* = { q L. +1) ^m(V . n+1 ) > 1 j-g 7 ^ 

\ m(V) + l if m(V 0;n+1 ) = 



To avoid confusion, indicate with the symbol ajj the Kronecker delta. Finally, abbreviate = (gj" 



air Willi I lie s\ n ii i( ii !(; , I i it' i:\ i ( uu tkcj iirua. riiianv. a i a )i ( '\ la ic q = \q 

(n+1) -, „(n+l) 
The bulk of the theorem is contained in the following assertion. 



Proposition 2. For any T> £ A(x n+1 ; [0,£]),1 < n < N — l,t > Q anrf almost all x n S r„, t/ie integral 
I(V)(x nl t) is given by the following sum of values of trees of A(x„; [0, t]): 

W = «m(B 0i „ +1 ),o {N-n- m(V)) V(V /0 . n+1 ) 

l* n+k-l 

+ E E V(v /0 . n+1 o k j +1) .v ., n+1 ) . (5.8) 

fc>l i>l 

Using the terminology introduced in Section TS Al we can give the following graphical picture of Proposition 
[2J The nodes divide each line of a tree in segments that we shall call branch intervals. To compute 
I(V),V=(x n+1 ,S)GA(x n+1 ;[0,t]): 
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P N (K N ,t) = 

—N 



P*-,^'*- f dx » £■ 




j=1 





N-3 /*5 N-3 // N-3 N-3 / N-3 // N-3 N-3 N-3 _ 

> 2 ' 1 >2' 1 >2' 1 ! 2' 1 >2 m1 >2 m1 4"' 

>3' 1 

etc. 



FIG. 7: Integration of degrees of freedom: from Liouville equation to BBGKY hierarchy. 
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1. extract from T> the subtree Po;n+li what is left is C/o ; n+i- 

2. Reattach 2^o;n+i through its root to any branch interval of 2?/o ; n+ii with the following care. The 
resulting tree will have the old m{T>) nodes of the starting tree T>, plus one new node to which the 
root of 2?o ; n+i is attached: the reattachment must be done in such a way that the reciprocal order of 
the old nodes in the resulting tree is the same as it was in the original tree T>. 

3. If the new node lies on the root line, append to it a node label with value in the set (1, • • • , n). 

4. Sum all the possible resulting trees found in points 2,3. 

5. If m(2?o;n+i) = 0, add to the result of point 4 the tree obtained by discarding 2?o ;n +i = T, i.e. 
D /o-.n+i = (x n , 5), multiplied by a factor (N — n — m(T>)). 

Several examples are provided by Figure [7] 



2. Proof of Proposition^ first step 

In this section and in the following we give an outline of the proof of Proposition [5J we refer to the 
appendices for the details. 

We shall begin with a couple of general definitions that will be useful along the proof. Let T> £ 
A(x n )[0,t]),T> = (x n ,5), with S as in (|3.3|) . Let k £ (1, • • • , m(T>)), and denote "D/w k the collection of 
variables T> deprived of variable w k . Similarly, (£x>)/fc(s) ((A/x>) / k ( s )) will indicate the state (number) of 
particles of the evolution associated to T> at time s, forgetting (when it is present) particle n + k. We define 
two subsets of the unit sphere surface by 

^*±(£>/tf>J := |ui fc such that V £ A{x n : [0,i]) and 

((£v)/k(h±s),T^(qj k (tk}T>)+aw k ,p k f) £ T { ^ v)/kitk±s)+1 Vs e (0,r ± )| , (5.9) 

where r+ = t — tk and r_ = t k — t c»+*) (the condition of existence of the free dynamics is understood in the 

ii 

definition). That is, Wk £ ^ < fj±('^ > /w k ) ensures that particle n + k moves freely with no collisions with the 
other particles of £t> for all times (up to t) in the future (+ case) or up to the time, in the past, in which a 
new particle is created through a collision with n + k (— case). 

For future convenience, given (x n ,5) £ A(x„; [0,f]), we introduce the set of collision histories for fixed 
tree 5, Ag(x n ; [0, t]), which is in one by one correspondence with the set of times, momenta and unit vectors 
with the constraints explained in the definition of collision history. We use the notation 

S = (ti, . . . ,t m ,pi, . . . ,p m ,W!, . . . ,w m ) , (5.10) 

Ajfen! [0,i]) 3 V = (x n ,5, S), and we define a measure dS on Aj(x n ; [0,i]) as the Lebesgue measure with 
respect to the variables ti, . . . , t m ,pi, . . . ,p m ,wi, . . . , w m . The value of the tree can be written 



V(V)= / ddR(x nl S,S). (5.11) 

"'A J (x ii ;[0,t]) 

The proof consists of two steps. In the first one, by using a careful subdivision of the integration region as 
well as appropriate changes of variables, we derive a formula which is the same as (|5.8[) except for the fact 
that the integrals over collision histories in the right hand side (hidden in the definition of V) are restricted 
to certain subsets. Namely, in the assumptions of the Proposition and using the notations 

T> = (£„+i , $) , 

V/o-n+l ° M ("+i> ;l £>0;n+l = (Xn,1k,i) , (5-12) 

we can prove 

I(P) = a m CD ., n+l) ,o (N-n- m(V)) V(V /0 , n+1 ) 

I* n-\-k — l „ 

+ Y1 d%iR(x n ,y k>i ,%i) , (5.13) 

k>i i>i J ^ k \(z n ;[a,t]) 
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for almost all x n G r„, where 

^iMn, [0.*]) ~ U 4t-(£n! tM) ( dis j° int Uni0n ) 

4t±^n; [0.*]) - {ft.i G A^.ifen! [0,*]), with 

such that w k € ^*( 7M ) ± (fe„,7fe,i>7fe,i)/iDj} • (5-14) 

In what follows we shall explain briefly and without formulas how this equation is derived, restricting for 
simplicity to the case m(T> Q . n+ i) = 0; the rigorous proof is in Appendix [B] In I(V) we integrate over x n+ i 
an expression which is given by an integral over the collision histories T> that are compatible with the tree 
T>, of a function R which depends only on the evolution S-p - see (|5.5|) . (|5 . 1 1 [) . After interchanging these 
integrations, we can parametrize x n+ \, that is the state of particle n + 1 of the evolution at time t, with the 
state of the same particle outgoing its last (in [0, i]) collision in S-p, when this collision exists. Such a state 
is described by the time of the last collision t* , the index i indicating which particle undergoes the collision 
with the n + 1— th, the unit vector w* := a~ 1 (q„ +1 (t*; V) — qi(t*;V)), and the momentum p* of particle 
n + 1 outgoing the collision (which is equal to p n +i)- Then, for x n +i such that this collision exists, say with 
i, we change variable — > (t*,p*, it)*): the resulting integrals J dt* J dp* J <») dw* correspond to a new 

node that has to be added to T>, while the Jacobian determinant produces the associated weight factor. The 
net effect is the value of a tree produced via operations 1 and 2 of the list at page 1181 when the integrations 
associated to the new node are restricted to "outcoming collisions producing a particle that does not collide 
with the particles of the evolution for all times in the future up to time t" . 

We are left with the integral over x Jl+1 such that the last collision does not exist. There the integrand is 
composed by a weight function which is independent on x n+ i (since we are assuming also m(V . n+ i) = 0), 
and a time-zero correlation function P n+1+m (x5) which depends on x n +i only through its correspondent value 

at time zero T_^ + (x n +i)- Therefore, by changing variable x n +i — > x 'n+i = ^-t+( x n+i) an d extending the 
integration over the whole one particle phase space compatible with the state of the other particles of £z>(0), 
we eliminate completely the particle n + 1 and recover a correlation function p n+m ^ , multiplied by a factor 

N — n — m(V) (see definition (|2.20p ). This correspond to operation 5 of the list at pagefT8l and produces 
the term in the first line of (I5.13[) . 

The error term in the preceding extension of the integration region will contain an integral over the states 
x' n+1 of the particle n + 1 at time zero such that "there exists a first collision (in [0, i\) with at least one 
of the particles of the evolutions associated to X>o;n+i" • This integral function is in turn integrated over 
all such evolutions, that is over all the collision histories V that are compatible with the tree Z?o;n+i- 
Calling t*,p*, w* the time, momentum and unit vector variables describing, in the usual way, the state of 
the free evolution of x' n+1 ingoing its first collision with the particles of Bqj'i we can proceed as before by 
making the change of variables x' n+1 — > (t* ,p* , w*). Again we have resulting integrals J dt* J dp* dw* 

corresponding to a new node to be added to £>o;n+i, and a Jacobian determinant producing the associated 
weight factor. The net effect is the value of a tree produced via operations 1 and 2 of the list at page IT51 
when the integrations associated to the new node are restricted to "incoming collisions producing a particle 
that does not collide with the particles of the evolution for all times in the past from t* up to 0" . This term, 
together with the one obtained in the above paragraph, gives, once summed over all the possible choices of 
particle i, the term in the second line of (|5.13[) . 

The case m(Z?o;n+i) ^ is treated in the same way, with the only important difference that the role played 
by time is now played by ^ = ^ , see Appendix lB 21 Here we only mention that, in particular, the 

"last collision" has to be understood in the time interval [ti*,t], and in the terms with "no last collision" we 

perform a change of variable x n +\ — > x' n+1 = T^ +t m (x n+ \). After this change of variables, the extension of 
the integration region to the whole one particle phase space compatible with the state of the other particles 
of £x>{ti*) gives a term that is shown to be identically null, using cancellations between outcoming-incoming 
collisions occurring at time ti*. This explains the Kronecker delta in the first line of (|5.13[) . 
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3. Second step: cancellations between collision histories 



Let us come now to the second step of the proof of Proposition [2j we will show that we can extend to 
^7 fc iilln'j [0, *]) the integral in the right hand side of (|5.13[) . since the total contribution of the missing set is 
equal to zero. This is achieved by the following 



for almost all x n G r„ . 

Formula (I5.15|) . together with (|5.13l) and (|5.11|) . complete the proof of Proposition^ 
We give the proof of the lemma in Appendix [C] and outline it briefly in the following. Consider the set 
of all the possible trees of the form Qk,% — (x n ,j k i ) — P/o;n+i ° fc 9 < ll + 1 ' i ^0;n+i, obtained from a given V 

through the operations 1 - 4 in the list at page [15] (A; is the node of Gk,i that is created in such operations, 
and jk(Qk,i) = *)■ In the left hand side of (|5.15[) we sum and integrate over the corresponding set of collision 
histories Qk,i such that the associated evolutions satisfy one of the two following special recollision properties: 

A. the particle generated in the node number k is in outgoing collision at time tk with particle i and, if 
we let this particle evolve forward in time together with the particles of the evolution £g k i (s),s > tk, 
it undergoes a new collision within time t; 

B. the particle generated in the node number k is in incoming collision at time tk with particle i and, in 
the backwards evolution £g k i (s), s < tk (of which the particle takes part), it undergoes a new collision 
within the time of the first node that is found climbing the line generated in node fc, or within time 
if there is no such a node. 

We shall refer to the "new collision" as the "recollision" with the other particles of the evolution. 

The lemma follows from the observation that for any evolution (i.e. collision history) of type A there is an 
evolution (collision history) of type B (and viceversa, so that a one by one correspondence is established), 
giving opposite contribution to the left hand side of (|5.15p . To find it, add to the evolution £g k t of type A 
the free flow of particle n + k from the time tk up to the time of the recollision (or, if you start from an 
evolution of type B, erase it from the time of the recollision up to the time tk)- Clearly the new evolution 
(and associated collision history) that is obtained in this way, say £g , .,, is of type B (or A), with k and i 
substituted by some different values k' < k(> k), 1 < i' < n + k' — 1. The two evolutions will correspond, in 
general, to different trees. Moreover, the value of function R(G k >,i>) is obtained from R(Gk,i) by substitution 
of the weight factor of node k of Qk.% with the weight factor of node k 1 of Qy ,v '■ But this is, up to a minus 
sign, the transformation in the integrand function induced by the change of variables (tk,Wk) — > (tk',Wk'), 
where Wk,Wk' are the unit vectors labelling the nodes fc and k' in the two collision histories. Hence the 
lemma is proved performing this change of variables in the restriction of the integral to the collision histories 
of type A (or B). 



Lemma 1. In the assumptions of Proposition^ and with the notations of i5.12\) , \5. 1$ , it is 




l* n+k-l 



d% i R(x n ,'y kji ,% i ) = 



(5.15) 



C. 



Sum over trees 



Using Proposition [2l Equation 



(|5.4[) becomes 




1 





N — nl 



E'eAfe^ilO,*]) 



I* n+k-l 



+ 





(5.16) 



veA(x 



x n+1 ;[0,t])k>l i>l 



for 1 < n < N, a.e. in P 
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The content of the square brackets is graphically represented by a sum of trees of A(x n ; [0,t]) with certain 
multiplicity factors. Hence, to deduce the assertion of Theorem [TJ we are left with the problem of showing 
that in this sum we have exactly N—n copies of each tree of A(x n ; [0, t}) - see Figure [7] for some example. This 
follows easily from analysis of the extractionfe growth operations described by Proposition [51 as explained in 
what follows. 

Fix Q € A(x n ; [0, with number of nodes < m(Q) < N — n. If m(Q) < N — n, consider the tree 
V — V [xn+x) G A(x n+1 ; [0,t]) which is obtained simply adding a coordinate x n +\ G Ti(x n ) (and such 
that x n+1 G r* +1 ) to the root of Q. It is m(2?Q°^ +1 ) = 0, and we see that J Fi ( x j dx n +±T^ a \x n +i) produces 

all the N — n — m(Q) copies of Q that can be obtained through operation 5 in the list at page [TBI 

Now we want to find all the trees in A(x n+1 ; [0, t\) that produce one or more copies of Q via operations 1 
- 4 of the list. Suppose that node number k of Q is created in point 2. Then it is clear that we have one and 
only one tree producing Q: this tree can be reconstructed with the following operations: 

1'. prune the subtree Gk and delete the node together with the label (if any) attached to it; 

2'. reattach the same subtree to G /k in such a way that the root line of Gk is superposed to the root line 
of G /k, and the reciprocal order of the nodes in the resulting tree is the same as it was in G\ 

3'. add a coordinate x n+ i G ^i{x n ) to the root, as well as a label n+ 1 to the new nodes crossed by the 
root line. 

We shall call V = V {x n+ i) the result of these operations. By construction, we have V G 
A(x„ +1 ; [0, t\) (for x n+l G r* +1 ), and § r , x ^dx n +\D {x n +\) produces a copy of G when node number 
k is created in operation 2 of the list at pageTl8l 

The operations 1' - 3' can be repeated for any node of G-, giving m(G) trees (some of which are possibly 
equivalent) in A(x n+1 ; [0, In particular, we have exactly m(G) different ways to produce G through 
operations 1-3, hence operation 4 gives m(G) copies of G- These copies, together with the previous 
N — n — m(G) copies obtained by operation 5, give the total number of N — n copies of G appearing in the 
square brackets in (|5.16[) , thus concluding the proof of the Theorem. □ 



6. CONCLUSIONS 

In this work we discussed a derivation of the series expansion used by Lanford Q to perform the 
Boltzmann-Grad limit, expressing the time-evolved n— points correlation function in terms of the higher 
order correlation functions at time zero for a system of N hard spheres in a finite volume. We established 
a new method of construction of the series based on step by step direct integration of degrees of freedom 
from the solution of Liouville equation, rather than iteration of the BBGKY equations. Each term of the 
expansion was written in the form of integral over some fictitious evolutions of particles called "collision 
histories" , for which we could introduce a convenient graphical representation. We showed that these graphs 
can be used to control the integration procedure leading from the expansion for p n +i to the expansion for 
p n . Mutual cancellations between collision histories showing special "recollision properties" were exhibited 
as an important part of the proof. The method provides a construction of the series expansion in a fixed 
full measure subset of the phase space, under the only hypothesis of some integrable bound for the density 
of the initial measure, and symmetry in the particle labels. This strengthens the results previously obtained 
in literature. We stated also an extension of the main theorem to initial measures with non definite number 
of particles. 

Without assuming continuity along trajectories of the initial measure, we could resum the final expansion 
and recover the usual BBGKY hierarchy of integro-differential equations for hard spheres, as originally 
deduced by Cercignani in Q (and rigorously obtained in Q by using the continuity assumptions). In fact, 
the final expansion can be also seen as the series solution of the Cauchy problem for the BBGKY hierarchy: 
actually this is the way it was presented in @ where, nevertheless, a rigourous discussion was still missing. 
In the hard sphere systems, the rigorous analysis of the dynamics and derivation of the BBGKY equations 
is more complicated than for smooth potentials (which was well known at the time of Q), because of the 
singular character of the interaction. Such an analysis was realized, for the Hamiltonian dynamics, in 2] 
and [9(, while the rigorous derivation of the hierarchy was first made by Spohn in [12J. In what follows we 
make some comparison with that note. 
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The starting point in [12] is an equation (Proposition 1, formula (20)) expressing the variation 

p n (xi,-- - ,x n ,t) — p n {T_t (x\, ■ ■ ■ ,x n )) as a sum over the number of collisions in [0,i] between the clus- 
ter of particles (1, • ■ • , n) and the others (n + 1, • • ■ ,N), of the corresponding gain and loss terms expressed 
through the initial probability measure P. This equation can be considered already as a rough form of the 
BBGKY hierarchy: the goal is to show that the sum of all the gain and loss terms is absolutely continuous 
with respect to the Lebesgue measure and has a density given by the collision operator of the hierarchy 
applied to p n +i, i.e. to compute the derivative of those terms with respect to the time. To do this, some 
probability estimate on the number of collisions is needed, together with the continuity along trajectories of 
the initial measure. After that, the expansion of Lanford is derived, as usual, via iteration of the hierarchy. 
Finally, it is rephrased as an integral over collision histories and then extended to non continuous measures 
of grand canonical type by a density argument. 

Coming back to the starting point, formula (20), we see that it does not keep track of what the external 
particles colliding with (1, • • • , n) do during [0,t], and that a control on the number of collisions in the time 
interval is required. As we saw, in the notion of collision history every time, going backwards, an external 
particle collides with (1, • • ■ , n), we add it to the cluster (1, • • • ,n) and keep looking at it. We saw also that, 
using this notion, we can directly express the variation p n {x\, ■ ■ ■ ,x n ,t) — p n {T_^ {x\, ■ ■ ■ ,x n )) as a sum 
over the number of new particles that can appear in the history, rather than over the number of collisions. In 
this way, provided we introduce the notion of collision history from the beginning, we can construct directly 
the final expansion (without the need of strong estimates nor continuity assumptions). Notice also that 
this construction is carried on through nothing but the same kind of steps leading to Eq. (20) of [T3|: a 
decomposition of the phase space, a flow of the coordinates (change of variables) from time t to the time of 
the last collision or to time zero, cancellations between sets. 

There are various other rigorous discussions on the hierarchy and the series expansion for hard spheres. 
A derivation of the BBGKY hierarchy in the form of integro-differential equations is given in [6| . There the 
authors show that, using the special flow representation introduced in [9(, a weak version of the hierarchy 
(integrated against test functions in a suitable space) can be derived. The final result follows then from 
uniqueness of the solution of the weak equations in the case of initial measure continuous along trajectories. 
Another discussion that has to be mentioned is in the work by K. Uchiyama [ijj], where the same results 
of [l2j are proved in a similar way, with the continuity along trajectories of the initial data substituted by 
the continuity over almost every point of the phase space (at the end removed again by density). Finally, 
other rigorous analysis on the Cauchy problem for the BBGKY hierarchy can be found in the works of D. 
Ya. Petrina and V. I. Gerasimenko, see Q, [l(| or the book [TT| . 

We hope that the methods presented in this paper can be used to deal with different situations. For 
instance, we believe that the whole analysis can be applied to discrete initial measures. Another direction of 
research would be the derivation of the smooth potentials case: following the ideas of [3] the procedure valid 
for the hard sphere case can be probably extended. Finally, it would be interesting to apply our methods to 
the construction of the series expansions in cases with boundary conditions different from those used here, 
and suitable for modeling of open systems. 
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Appendix A: Dynamics of hard spheres 

In this appendix we state the properties of the dynamics of hard spheres that, together with Proposition 
Q]of Section [2 Ai are used in our discussions. 

First, we prove formula (|2.15|) . We recall that the set n K n is the maximal subset of the n— particle 
phase space over which our main result can be derived pointwise for all times. We do not know whether 
r* = rt . However, we have 

Lemma 2. For any n < N, the set T n \ has Lebesgue measure zero. Moreover, dT n n (T n \ T' n ) is null 
with respect to the induced measure over dT n . 

The induced measure over the boundary da is given by Eq. (|2.8[) . The lemma is a simple consequence of the 
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existence of the dynamics over the full measure set T* , stated in Proposition!]] Since we have no information 
on the structure of r„ nor its measurability properties, we will prove the lemma via abstract arguments. 

Proof. It is sufficient to prove the assertion for any finite bound on the energy. A little abuse of notation 
will be used in what follows: we indicate with the usual symbols r n ,r* ,T n ... the bounded sets corresponding 
to an energy of the system not larger than E > 0, and with • |, | • |* respectively the restriction to the various 
sets of the Lebesgue measure and of the usual Lebesgue outer measure (that is the infimum, over all the 
possible coverings of a set built up with n— dimensional boxes, of the sum of the measures of such boxes). 
In the following we will use that the flow of the dynamics preserves also the outer measure, which can be 
easily deduced from the fact that it is an invertible measure preserving transformation (see [9j, p. 651). 

We abbreviate 

z (0) r * \ pf(O) 

z„ = r*\rt = |jTi")(4 )). (A.i) 

By Proposition [U it is |r„ + fe \ r* +fc | = for any k, hence by the very definition Z n ^ must be a null set. By 
the same reason, it suffices to prove that Z n is a null set too. To do so, we use a contradiction argument: 
suppose that Z n is not null; we will show that this implies the existence of a not null subset of T n+ k \ ^ n +k 
for some k > (which is forbidden by Proposition [lj . 
For any x n G T* , 1 < k < N — n, call 

B k {x n ) = {y k G T k (x n ) s.t. (x n ,y k ) 6 r„ +fe \ T* n+k }, (A.2) 



so that we can write Z n ^ = UfcLi™ %nk> 

Z [ °i = {x n G r* s.t. \B k (x n )\* > 0} . (A.3) 
Given a function rj : zh°^ — > (0, oo), define also 

N-n 

U U ^){x n ), ^=\JZ n % (A.4) 

Observe that, in the assumption \Z n \* ^ 0, there exists necessarily a value of k such that, 

for any V (x n ) > 0, \Z^ k \* > . (A.5) 

Otherwise, take 770(2,1) > for which \Z n Vo k \* = for all k, and let e m be a sequence of positive 

numbers converging to zero: writing Z n = \Jm=i Ujez^i^xJ € and ^ofe„) > £ m } = 

\JZ=i\Jj& T $}L ( Z n m) r\{Tl n \x n ), with x„ G Z£>\r) (x n ) > s m and s G K}) , we would get \Z n \* = 

by subadditivity and preservation of the outer measure (the set in the argument of T^jg is a subset of Z n , 
hence it has outer measure zero). 

From now on k indicates the variable for which the condition (|A.5|) holds. Given a function 77, we can 
consider the following subsets: 

B k (x n ) = {y k G B k (x n ) I 3 T^ +k \x n ,y k ) Vs G [-jjGO.t/GeJ] 
and T^ k \x nl y k ) = (rWfe,),T(\))} , 

r<f>) -Itn. ..\*+ ~ a r.W 



= {(x n ,y k ) s.t. x n G ZZi and y_ k G B k (x n )} , 
W%1= [J U T^(x n ,y_ k ) 

= lj U (T^(x n ),T( k \B k (x n )j) . (A.6) 



£Z (o; 1 s g[-,( £ J,,(jJ] 
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By definition W^ k C r„ +fe \ r; +fc , and since zi° ] is null we have \W^ k \ = 0. Points of W^ k do not have 
well defined evolution for all times, but still the evolution exists up to times r](x n ), and this enables to define 
W^l . Notice now that for any x n G Z^ 1 

\Bk{x n )\* > \Bk(x n )\* — \{y k 6 Fk(x n ) such that they experience 

at least one collision between themselves or with particles in x n 

or with the walls, when evolved with the (n + k)— th particle dynamics, 

within time [-v(x n ),v(x n )}}\* 

> \B k (x n )\* - 0( V (x n )) > \ Bk{ £n)\* > {A 7) 

for rj{x n ) sufficiently small (the bound with 0(j)(x n )) can be obtained by simple geometrical estimate; see 
for instance Q, p. 24-26, which can be easily adapted to our case). 

Choose a function rj(x n ) such that the above inequality holds. Then from the last line of Eq. (|A.6|) . using 
(|A.5p . (IA.7I) and preservation of outer measure, we see that it must be 

\W%1\*>0. (A.8) 

Since C T n+ k \ T* +fe , the contradiction is found. This proves that r„ \ r£ is a null set. 

To prove the second assertion of the lemma, notice that a not er„— null set A n over dT n in which the 
dynamics is everywhere well defined, spans a set of strictly posive outer measure over r„ through the 
operation U sg [ T ],T > 0. In fact, the time return to <9r„ is r(x n ) > for almost all x n of the set, and 
the Lebesgue measure over the subset of points of r„ whose previous collision was in dT n is da n dt, t being 
the time elapsed after the collision ([§], [j]). Hence each "box" B n of dT n spans at least a set of measure 
J B da n (x n )T(x n ) in r„. Conversely, each box of positive measure in T n corresponds to a set of positive 
measure over dT n : we refer to @ for more details (see Lemma 3.1). Since lj se [ T AZ n C\dT n ) is a subset of 
Z„, which has been shown to be null, it follows that the set Z n n dT n must be also null in the measure da n 
over <9r ra . This, together with Proposition [TJ completes the proof. □ 

Let us turn now our attention to the set K. n . It is unclear whether JC n coincides with T^. In any case, 
what is relevant for our purposes is formula (|2.16|) . To deduce it, we state another known feature of the hard 
sphere dynamics, which is related to the collision surfaces. Denote with dX the Lebesgue measure on R. 

Lemma 3. Given a set A C T n with \A\ = 0, then T^ n \x n ) A for almost all (x n ,t) £ dT n x M, with 
respect to the product measure do~„ x dX. 

The lemma can be easily deduced from the properties of the special flow representation discussed in Q , Q • 
For a complete proof, we refer to [3] (Lemma 3.4). From Lemma [3] it follows 

Lemma 4. For any n < N, the set T n \ K, n has Lebesgue measure zero. 

□ 



Appendix B: Proof of (pTT3l) 

1. Case m(X>o;n+i) = 

Consider a tree V G A(x n+1 ; [0, t]), V = {x n+1 ,S) and suppose m(2? ;n+i) = 0. We will prove the 
statement in this case first. 

From (|5.5p and (|5. 1 1|) we see that 

I(P)= [ dx n+1 [ d5R(x n+1 ,5,S) 

JTtixJ JA T (x n + 1 ;[0,t]) 

= [ d5'R(x n ,x n+1 ,S,S) , (B.l) 

" / A„ + 1;I ( ;£ „;[0,t]) 

where we called A n+1 .g(x n ; [0,t]) the set of collision histories ^-> Xn+1 eTi(x )^j(£ n +ii [0, i]) (which is in one 
by one correspondence with the elements 5' — (x n+ i, 5)), and dd' = dx n +\d5 the measure over this set. By 



27 



assumption (|2.17l) - see also (|3.12p and discussion above -, R is a summable function over A n+1 .g(x n ; [0, t]), 
and all the integrals can be interchanged freely. 
Now introduce the subsets 

A i+i ; lfe™; [0,*]) := {V G A n+1 .j(x n ; [0,t]) such that ^ n+1 (s;P) = rl^+^+i) Vs G (t,.^,*)} , 

^Jfe- t '^) := ^ 6 A n+i;5fe„; [0,t]) such that x„+i( S ;P) = TlV+.+ ^+i) Vs G (t*,t) , 
** G (tfc,*fc-l) and g n+ i(i*;X>) - tfc(**;D) = aw*, 

\w*\ =l,w* ■ (p n+1 (t*;V)- Pi (t*;V)) > 0} , (B.2) 

for 1 < fc < l*(S),l < i < n + k — 1, where Z*(<$) is the variable defined in (|5.7I) . In our assumption 
= m(5) + 1 (and t t ,^ =0): we give the definition in this way because it will be useful to deal also with 
the more general cases. We remind the reader that the configuration of a particle in the evolution associated 
to a collision history is defined as the limit from the future of the flow of the dynamics; for example in (|B.2|) 

it is q n+ i(t*;V) = T ( _^ +t , + (x n+ i), etc. Then 



l(V) = / dS'R(x n ,x n+1 ,S,S) 

',*]) 

/ dS'R(x n ,x n+1 ,Sj) . (B.3) 



i(<5) + l n+fc-1 

E E 

k=l l= l jA ZW^nA0,t]) 



Put, as usual, D = (x n+1 , 5, 5). In each term of the sums in the second line of (|B.3|) we can perform the 
change of variables 

x n+1 ^(t*,p*,w*), (B.4) 

where t*,w* are the variables introduced in the definition of the integration sets (|B.2|) . and p* := 
p n+ i(t*]V) = Pn+i- That is, t* is the first time of collision of particle n + 1 with the other parti- 
cles of the collision history going backwards in time, particle i is the one colliding with n + 1, and 
w* := a -1 (q n +i(t*; T>) — qi(t*;T>)). Then it is a simple exercise to see that the measure transforms as 
dS' — a 2 w* ■ (p* — Pi{t*;V))dt*dp*dw*d8 (see for instance the Appendix 4.B of Q), where pi(-) does not 
depend on the full T> = (x n+1 , <5, 5) but just on (t*,x n , S, S) = (t*, £>/o ; n+i7 S). Rename the dummy variables 
as (ti,pi,wi) — >■ (U + i,pi + i,wi +1 ) for I = k, k + l,--- ,m(S), and (t*,p*,w*) — > (tk,pk,Wk), and call the 
new resulting set of variables 

%i ■= (ti,-* - »**!,'•• »* m (5)+i.Pi'-- - >Pk,'-- »P m (I)+i>^i.--- >fe,-- - ,M> m( i)+i) , (B.5) 

and also := dt*dp*dw*dd. 
Consider now the tree defined by 

Qk,i = (£n,7fe,i) : = ^/o^+i 2?o ; n+i G A(x n ; [0,t]) (B.6) 

(in our case it is V / . n+1 = (x n ,8), and X>o;n+i = 7"), and consider the collection of variables 

Qk,i '■= fen'7fc,i>7fc,i) • (B.7) 

We shall see that, at least for a. a. i„ G T* , the domain of integration of the new variables is the set of jk t i 
such that Qk,i is a collision history in [0,t], with only one additional constraint on Wk, which implies that 
particle created in the outcoming collision at time t/. would move freely in the future (since in (|B.2[) t* is the 
first (backwards) time of collision of particle n + 1 with the others). 

First of all, it is clear that we can assign to Qk,i an evolution £g k i (s),s G [0, t] , in the same way as we do for 
collision histories, and that this evolution is well defined in our domain of integration for almost all x n G T n 
(the evolution coincides with £x>: just erase particle n+ 1 in the time interval (t*,t)). Then, our integration 
region is defined as the set of such that: (i) < t m( ^ .)= m (S)+i < t m(~ k < ■ ■ ■ < h < t ; (ii) 

Pi,--' >Pm(7 feii ) G r3 ; (iii) for 2 = ,fc-l,u)j such that ^n+i-i^S^i),^^^^;^^) + awi,pi) G 
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T n+ i ; (iv) Wk G fi^T {(Gk,i)/w k ) ; ( v ) Wk is such that the clusters of particles of Gk,i, (1, 2, ••■ ,n,n + 
/c), (1, 2, • - • ,n+l,n + k), ■■■{1, 2, ••• ,n + k-2,n + k) are respectively in T* +2 , •• • ,r* +fe _ 1; i.e. they 

do not run into singular configurations; (vi) for I = k + 1, • • • ,m( 1 J ki ),wi G ^(7*. &k,i),Pl) ■ 

Now, restricting to x„ G T*. consider the difference between the set defined by (i),...,(vi) and the set 
Al* , + {x n ; [0, i\) defined in (|5.14[) (and equal to such that conditions (i), (ii), (iv), (vi) hold, and 

condition (iii) is modified by replacing r„ + ; with that is by ibi G ^',(7. )(^n+i-i(*') Gk-,i),Pl) for 

I = 1, ••• ,k— 1.}); this difference contains only values of 7^ such that some subcluster of particles of 
(1, • • • , n + k) run at some time into a singular configuration (and we can also notice that this singular 
configuration does not occur, in any case, along £g k 4 (s) for s G [0, t\). Hence, for almost all x n G T* the 
integral in dr/k,i over the difference set must give zero contribution: otherwise we could find, in the phase 
space of such cluster of particles, a set with Lebesgue measure different from zero over which the dynamics 
is not well defined (contradiction with Q, [!])■ We do not give a formal proof of the last statement (which 
is not difficult to believe): this can be found in [l4[ (see Lemma 6.2 of that work). 

In conclusion, noticing that, after the above renaming of the variables, Pi(t*; f/ojn+ii = Pi(tkl Qk,i) and 

a 2 w* ■ (p* - Pl (t*;V /0 . n+1 J))R(x n ,x n+1 (t\p*,w*),6,6) — ► R(Q k>i ) , (B.8) 

we have obtained 

/, , dS'R{x n ,x n+1 ,d,S) = / d% ti R{x n ,i ki ,*( k<i ) (B.9) 

JA^ip fe„ ; [0,t]) J A^ . + (x n ;[0,t]) 

almost everywhere in T n . 

Now we want to deal with the term in the first line of (|B.3j) . For almost all x n G T* the integration region 
in the term considered can be rewritten as 

1 5' = (x n+1 ,5) such that (V /0 , n+1 ,S) G A ? (x„; [0,t]) and 

(^ / o i „ +1 ,a)( s )» T -t ) +s+ (^+i)) e r^ ( _ oib+ii , )W+1 v« g (o,t)} ; (b.io) 

in fact we can notice, as done just above, that the error term is an integral over a region of zero measure, 
corresponding to singular trajectories. In the region (|B.10[) the dependence of the integrand on the vari- 
able x n+ \ is concentrated on the correlation function, since the particles of the evolution appearing in the 
definition of the set evolve independently of x n+ \ in our assumption m(2?o;n+i) = 0. Explicitly, 

R(x n ,x n+1 ,6,S) = W (v /0 , l+1 ,d) Pn+l+m{1) (T%{x n+1 ),£^ fQin+u ~ 5) {Q)) (B.ll) 

(here we used the symmetry of the correlation functions). 
Hence by making the change of variables 

x n+ i — > x' n+1 = T { }} + {x n+ i) (B.12) 

we obtain the integration over 

|<5" = {x' n+1 ,8) such that (V /0 . n+1 ,S) G A 7 (x n ; [0,t]) and 

(^n+i^W^i- G T ^ /0 ., n+1 M°^ VS G ' (B.13) 

in dS", of the function 

W (p /0 . n+1 j) p n+1+m( - s) (<+i^ ( p /O; „ +1 i)(0)) ■ (B.14) 

We want to complete now the integral in order to obtain the function p n+m (gy This can be done extending 
the integration to the full set 

{S" = (x' n+1 ,S) such that (V /0 . n+1 J) G Aj(x n ; [0,t]) and x' n+1 G ^^(^^^(O))} . (B.15) 
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The integral in dS" , over the region (|B.15|) . of function (|B.14I) gives, after ordering the integrations, 
dSW (v /0 . n+1 ,S) / dx' n+lP { g- (x' n+1 ,£(p t s)(°)) 

v ' J rn+~o ) V(p /0sn+1 ,i ) m [ ' K /0 - +1 ' ' 

= (N-n-m(6)) f ^ dSW (v /0 , n+1 , s) p n+m( - s) (%> /Oin+li j)(0) 



A 7 (z„;[0,i]) 

= (JV - n - m(J))7(© /0in+1 ) • (B.16) 
Subtracting the error term, we have obtained 

dS'R{x n ,x n+1 ,5,8) = (jV-n- m(5))V{V /0 . n+1 ) (B.17) 



m(S) + l n+fc-1 

E E 



fc=l i=l ^i^fenSP.*]) 

for almost all x n € T n , where 



^(fc,t;- 



dl"W (p/0;n+l,<5) P n+1+m( 5) ( a! n+l> £ (5 /a . B+1) «)( ) 



n+1; 



.y^i [0,*]) := [S" = K +1) S) such that (x n ,S,5) € Aj(x n , [0,t]) and 

t* e (tk,tk-i), and (T a ^?(a4 +1 )) g - g» (**;£„, 6, 5) — aw* , 

|t&*| = l,t&* ■ ((T, W (x; +1 )) p -Pi(t*;x n ,S,6)) < O} . (B.18) 

In expression (IB.17[) we have decomposed the error term in a sum of integrals, where labels k and i describe 
between which nodes and with which particle of q +i ^ occurs the first collision, moving forward in time, 

of the external particle with initial configuration x' n+1 at time 0. Once again to write (|B.17[) we removed 
sets of zero measure (i.e. points x' n+1 6 Ti 

We can treat the terms in the second line of (|B.17[) as we did for those in (IB.3|) . In this case we perform 
a change of variable 



(t*,p*,w*), (B.19) 



where t*,p* are the variables introduced in the definition (IB.18[) and p* := p' n+ i- The measure transforms 
as dS" = —a 2 w* ■ (p* — pi(t*;x n ,S,S))dt*dp*dw*dS. We rename the dummy variables as (ti,pi,wi) — > 
(ti+i,pi+i,wi+i) for I — k, k + 1, ••■ ,m(8), and (t*,p*,w*) — > (tk,Pk,u>k), and we introduce the same 
notations of (|B.5j) . (|B.6j) . (|B.7|) . We can assign to Qk,i an evolution £g k:i (s),s <E [0,t], which is well defined 
in our domain of integration for all x n € T* , and that is obtained by adding to £(p +1 J) the free flow of 

x' n+1 in the time interval [0,tfc]). Moreover, for almost all x n € T*, the domain of integration of the new 
variables is the set of jk.i such that Gk,i is a collision history in [0,t], with only one additional constraint on 
Wk, which implies that particle created in the incoming collision at time moves freely in the past (since in 
(IB. 181) t* is the first (forward) time of collision of the particle starting in x' n+1 with one of the others); this 
is so by forgetting, as usual, the zero measure sets in which some cluster of particles of Qk,i run at some time 
into a singular configuration. This means that, making use of the definitions (|5.14[) , (|5.9[) . and rewriting the 
integrand with the notations introduced, 



dS"W {V/ a . n+ i,5) p n+1+m( - 5) (<+i>%> /o . n+1 ,,S)(0)) 

#Kifl(2 B ,7fc 1 i>7k.O (B.20) 



A^ 4 _(x n ;[0,4]) 

almost everywhere in T n . 

This last equation, together with (|BT7|) . (jB~3l) and (|BT9|) . gives ((57T5jl . 
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2. Case m{Vo-„+i) > 

Let us consider a tree V G A(x n+1 ; [0,t]),V — (x n+1 ,S), with m(X>o;n+i) > 0. This case is very similar 
to the previous one and it is discussed essentially in the same way, with the only difference that the role 
played by time t m ^ +1 = is now played by t>i*(g) = ^ q [ n+1) (S) ~ see Q5.6[) . (|5.7|) (through all this section I* 

and will indicate the values associated to 8 defined by (|5.6I) and (15.71) ). 

In particular, the analysis from (|B.1[) to (|B.9I) is exactly the same once we restrict to I* the sum over k in 
(|ET3j) . and substitute (|FTd| with 

= fen,7*,i) : = ^/0;n+i ° fe , £ <"+i) ;l ^o ; „+i £ Afe n ; [0,*]) . (B.21) 

Hence we have again 

I(V) = / dS'R(x n ,x n+ i,S,8) 

JA« ;I (2„;[0,il) 

+E E 

fc — 1 i — 1 Tfe i+ * — n ,L " ' J/ 

almost everywhere in T n . Of course now T> / . n+1 =^ (x n ,S), and 2?o:n+i 7^ 7", but they have a more 
complicated structure depending on the labels attached to the nodes on the root line of the tree V. 

Call 2?+' the collection of variables obtained from the collision history T> = (i n , x n +i, S, 8) by depriving it 
of x n +i and of the variables associated to the nodes with ordering number larger than I* — 1, and substituting 
m(6) with 1* - 1. With the usual notations (Eq. (133)) and (I5~TU|) ). 

= {x n ,l* - 1, Ji, • • • , ji*-xM, ■ ■ ■ - - ■ ,p;«_i,u>i, • • • . (B.23) 

Then, for almost all x n+1 G r n+ i such that 2? G Aj^-fjz;^; [0,t]), it is also 2?+' 1 G A(x n ; that is the 

same collision history restricted to the time interval and deprived of particle n + 1. For these values 

of x n+1 , putting x' n+1 = T_t +t t + (x n+ i) , we may define also a collision history with time span [0, U*] by 

V ( * ] = (£ v M(U*),x' n+11 q' n+1 +awi»,pi»,m(S) - I*, 

ji'+ir- J m (s)A*+i,-' - >* m (I)>^+i!-- - >P m (S),u>i'+i>"- ,ti m(1) ) . (B.24) 

This will belong to A^ M +2 (£ v (*> (U*), x' n+1 , q' n+1 + awi*,pi*; [0, *;*]). 

We can rewrite the integration region on the first line of (IB.22|) as the set of values of 
5' such that x n+1 (s;V) = T% s +(x n+1 ) Vs G (t,.,i), 25 J> G A(x n ; [t,.,t]) and V { * ] G 
Aaa , , +2 (£-.(«) (iz*) , i'ii, o' ■ 1 + awi*,pi*; [0,£;»]) (we are just discarding a zero measure set for almost 

all x n G r„). After that, we perform the change of variables 

Xn+i — > x' n+1 = TrJ +tl , + (x n+1 ) . (B.25) 

We obtain 

/ , , d6'R(x n ,x n+1 ,S,6) = [ dS" ( J] Wt(pP)) a 2 m> ■ (ft. -p' n+1 )R(V M ) , (B.26) 

where A is a short notation for 

A := {£" = « +1 ,<5) such that 25^ G A(z„; , 

G A ^ D (, )(t j+aC^w ^n+u + awi*,pi*; [0,fj.]) , 
and (^<„( S ),rW t+s _« +1 )) G r^ ( . )W+1 Vs G (t,.,t)} . (B.27) 



fa; ifO.tll 



d%,iR{x n ,~f ki ,%, 



(B.22) 
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Observe that the values of Wi* ,pi* in 6 associated to the particle colliding, at time ti* , with the one in x' n+1 

(*) 

in the evolution T>_ , describe both outgoing and ingoing collisions: we will strongly use this fact at the end 
of the proof. 

Extend now the integral to the integration region 

B := \s" = (x' n+1 J) such that V { * ] e A(x n ; [i,.,t]) , 

V { * ) eA AfM +2 {£ pW (t,.),< +1 ,«M-i + o«>i.,A.;[0,t,.])} , (B.28) 

and notice that the error term is an integral over the set of variables such that, for some l<k<l*,l<i< 
n + k — l,t* e (tk,tk-i), w* e S 1 , it occurs that 

(£ v m(s),T%+s-(<+i)) 6r v (>)+1 (B.29) 

- %(**;^) = «&* . 

i&* ■ ((Tl 1 t ) ! , +t »_(<+i))p -^(t*;^)) ^ (hence < 0) (B.30) 

(particle i is now identified ordering, in the usual way, the particles of V ( * ] G A(x n , [*/.,*])). 

Then in S\.A, callings* = (Tij^ ^ +t »_(x^ +1 )) p , we can associate to x' n+l the triple (t* ,p* ,w*). Adding 

to S the triple (t* ,p*,w*) and renaming the variables as explained before (|B.5[) . we obtain a collection 7^ 
defined as in (IB.5|) . This, together with (IB.7|) and (IB.21|) . defines a collision history associated to the tree 
Qk,i, as soon as the corresponding clusters of particles do not run into a singular configuration. Clearly, in 
this case it must be Q k>i e Ai^ ._(x n ; [0, *]). Moreover, dS" = -o?w k ■ (p k -pi(t k ;Q kA ))d-y kji (where t kl p k ,w k 

are now the elements in Gk.i)- By performing this change of variables and erasing sets of zero measure, we 
see that Eq. (IB.26|) becomes 



for all s S (ti* ,**), and 



l* -l 



[ dS'R(x n ,x n+1 ,S,S) = f dS" I J] ^(Df'lLv ■ (ft. - P ' n+1 )R{V { * ] ) 



IB 

r(s) n +k-i r 

^ Y] d%iR{x n ,*f ki ,%i) (B.31) 



fe=l 8=1 



for almost all i„ G T n . 

Furthermore, the first term in the right hand side of the above equation is equal to zero for almost every 
x n 6 T„. In fact, there exists an involution that associates to each element S'{ of B another element S 2 of 
the same set in such a way that the corresponding values of the integrand function have the same modulus 
and opposite sign. This involution is given by the collision rule applied to the two particles colliding at time 
ti* in T>_ , as explained in what follows. 

Given 6" = (yiji), yi = ((yi) q , (vi) P ), 

Si = - ,*/•,••• ,t m (s),Pi,--- ,P m (sy >«>{.,••■ ,W m (S)) > ( B - 32 ) 

by definition of 2?_ it follows that in its starting time £;» we have always a particle in the configuration 
((2/1)9 + au>j.,pV)- We put 5' 2 ' = (y 2 ,5 2 ) with (y 2 ) 9 = (yi) g , (V2) P = (yi)p + m*[wi* ■ (pi* - (yi) P )], and 
62 equal to Si except for the component pi* which is replaced by p\, — pi* — wi* [wi* ■ (pi* — (yi)p)]- The 
element 5 2 will belong to B. Looking at the integrand function, notice that the transformation S'{ — > S 2 
leaves unchanged the value of R(V^), as well as the value of the functions W r (V^) for all 1 < r < I* — 1, 
but transforms a 2 wi* • (pi* — (yi) p ) into a 2 wi* ■ (p[, — (y 2 )p) = —a wi* ■ (pi* — (yi) P )- Hence the integrand 
function changes its sign. 

Hence, equations (|B.22|) and (|B.31[) give the result. □ 
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Appendix C: Proof of Lemma [TJ 



Fix x n G r*, and look at elements %j G A- fc . (x n ; [0,t]) \ (x n ; [0,t]). By definition (|5.9[) and using 
the notations (IB. 51) and (|B.7[) . we have two cases: 

1. Wk G % fc (- fc i )+(£n+fc-i(^fc; Gk,i),Pk), and there exists 1 < k' < k and 1 < z' < n + fc' — 1 such that for 
some time t* G (ifc/,tfe/_i),u)* £ 5 2 , it is 

(7 fc ,i)(^' ^m) + aw k,Pkj^ G IV Sfc ,( s )+i VsG (£&,£*), 
( r -tl+t«-(9j fc (7 fc , s )(^;^fc,0 +a«)fe,Pfe)) -qi'(t*;g kt i) = aw* , 

((rS+t.-(*»(7 M )(**;^,0 + a**.^)) p -p<'(*^SM)) <°; (C- 1 ) 

2. Wfc G % fc (7 fc )-GE.n.+fc-i(*fej Gk,i),Pk), and there exists fc + 1 < k' < I* and l<z'<n + fc'— 1 such that 
for some time t* G {tw ,tk'-i),w* G S 2 , it is 

((^e fc ,J/fc( s )> r -t ) fe+s + (?i„Cr*, 4 )(*fc;^M) + o«>fc,Pfc)j e r^ Sfc s ( s ) Vs g (**,*&) , 
(r[ 1 l +t , + (q jk{ - k t) (t k -,g kti ) + aw kl p k )^ - qi>(t*;g kt i) = aw* , 



w 



Denote 1Z+ and 1Z- the sets of triples (k, i,7fc,i) with 1 < fc < Z*, 1 < i < rt + k — 1, and 7^ G 
A^ fc 4 fe n ; [0, t]) \ Ai*"* fe n ; [0, t]) satisfying respectively property 1 and property 2 of the list above. Introduce 
a measure dg over 7^+ U 1Z- as the counting measure with respect to k,i and the Lebesgue measure with 
respect to jk,i, and rewrite the left hand side of (|5 . 15[) in the short notation 



dQR(x n ,7k,i,%i)+ dgR(x n ,'y kti ,% t i) . (C.3) 

K + JTZ- 

We may define a transformation J over almost all 1Z+ as 

J[k,i,%i) = {k',i',T)k>,i>) , (C.4) 
where k',i' are defined in point 1 of the list above, and 

Vk',i' = (t'l, ■ ■ ■ >*m(7 fcti )'Pn ' ' ' , ^m(7 M )'"'l' ' ' ' '""'m(7 tei4 )) (C-5) 

is constructed from the collection of variables 7^, substituting the elements (tk,w k ) with (£*,?«*) (defined 
in point 1 of the list above), and reordering the components to obtain the usual decreasing sequence of times. 
Notice that (x„,7 fe / ^) is the tree that is obtained from j k i pruning the subtree generated in the fc— th node 
and reattaching it to the line representing particle i' , in such a way that a new node with ordering number k' 
is created, and that all the other nodes maintain the same mutual ordering. Then, it is clear that rjk',i' is in 
A- fc , ., (x„; [0, t}) as soon as all the clusters of particles associated to it do not run into singular configurations 
for all times. This is true for almost all jk,ii at least for almost every x n G T* . Moreover in this case, by 
construction, f)k>,i< is in A- fc , .,(x n ; [0, t]) \ aI, t (x_ n \ [0,t]), and it satisfies property 2 in the list above, i.e. 
(fc', i 1 , fjk',i>) is in 1Z-. Hence 

J : 1Z* + — > K- (C.6) 

where 1Z* + C 1Z+ and 1Z + \TZ*l has zero measure (for almost all x n G T*). Furthermore, the inverse function 
J -1 is defined over almost all 1Z- in a natural way. In particular, 1Z*_ = J(1Z* + ) C 1Z-, 1Z- \ 7Z*_ being a 
zero measure subset. 
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After substituting 1Z+ and 1Z- in (IC3|) with 1Z*, and TVL, we perform the change of variables (k, i, % t i) — > 
J(k,i,"jk,i) in the first integral. The function 

m (ik,i) 

Pn+m(-f k i ) fen+m(7 fc <) (Sk,i)) (C7) 

r=l 

is invariant under this transformation, while the measure transforms as 

d% tl W k (x n ,j ki ,j k:i ) = -d%',i'Wfc'fe n ,7 fc , )j /,d^ fe /,i') ■ (C.8) 
Therefore, the two terms in formula (|C.3[) cancel each other. □ 

Appendix D: Continuity properties 

In this appendix we prove some property needed in the discussion of Section 14 Al We always assume to 
work with an initial measure P with density /jv £ Cn', Liouville equation and correlation functions are 
defined by (|4.5[) . which is assumed to hold, for simplicity, on the whole set vl} + ^ (defined as in (|4.4|l ). The 
value of trees is defined by (|3.10jl . 

The following mild continuity property of the correlation functions can be derived with no need of addi- 
tional assumptions on the initial measure (and, as expected, V(T>) inherits the same property as a function 
of (x n ,t)). 

Lemma 5. For all x n e t}i + \ the functions of time 

t^ Pn (T t {n \x n ),t) , 

t^V(V)(T t (n \x n ),t) (D.I) 
with T> S A(a; n ; [0, t]), are continuous for all t > 0, that is 

lim p n (t}™ } (x n ),t + s)= lim p n (t!™ } (x n ) , t - e) , 

lim V(p)(T<g(x n ),t + e) = lim V(V)(T^(x n ), t - e) (D.2) 

e->0+ e— >0+ 

hold for all t > and all x n g . In particular, Eq. IIP. Sty is true for all t > and almost all x n G dT n , 
with respect to the measure do~ n . 

Remark. The continuity property stated in the above lemma is a consequence of the Liouville Equation, 
and it does not imply the stronger "continuity along trajectories", i.e. properties (i) and (ii) in the Remark 
(1) of Section QJ which are in general not valid unless we assume Eq. (|4. 3[) for the initial measure. 

Proof. We will deal first with correlation functions. For n = N the claim is a trivial consequence of the 
Liouville equation (|4.5j) . since the considered function is constant in time. 

Suppose that the property holds for the function p n +\ for some n < N — 1. Take e > small, and define 
T^ e \x n ) as the one-particle configurations x n +i compatible with x n and such that the evolution yj™ +1 ) 
docs not lead to a collision of the (n + 1)— th particle with the others in the time interval s £ (0, e]. Then 
we have 



dx n+1 p n+ i (T^l + (x n ),x n+1 ,t +ej - / dx n+1 p n+1 ( 
ri(T t ( ;» + (xj) V ' ^(TWfe,)) v ' 

Pn+l 



dx n+1 

r £) (T<;»( £ j) 



+0(e) - / dx n+1 p n+ i [T t ^\x n ),x n+ i,t) . (D.3) 

The term 0(e) is the restriction of the integral in the first term of the first line to the points that, evolved 
backwards in time together with the configuration of particles T^ s ,(x n ), display a collision with one of 
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these particles in the time interval (f, t + e]: explicitly it can be written, with the usual change of variables, 



a.s 



11 rS f n 

V" / dtx dpi / dw^wx ■ (pi -Pj(t + ti)) 

•p„ +1 (r^ + ( £ J,T a t ) 1+£+ fe(i + ii)+a*i,Pi),i + £) , (D.4) 

Pj (t + h) = (T t % + (x n )) , qj (t + tO = (T t % + (x n )) , and Q { + £ \- • • ) denotes the subset of 

Qj + selecting particles that do not collide with the others when evolved forward in times of the interval 
(f + fi, t + e]. Clearly the term in (|D.4[) goes to zero as e for e — > in our assumptions. 
A term similar to (ID.4[) is given by 



where here 



ri(T^j)\r<+ £ \T t ( ;>( £ j) 



dx n+lPn+1 (T^ +1) (T^\x n ),x n+1 ),t + e) = 0(e) . (D.5) 
Hence (|D.3I) becomes 

p n +i (r^^^Xn+O.t + e) -pn+i (T t ( | l) (x„),x„ +1 ,t)] I + 0(e) . (D.6) 



dx 



n+l 



Dominated convergence and the inductive assumption imply that this flows to zero with e. A similar analysis 
can be performed for negative e, therefore we have shown that, for all x n gT', 

fen p„(T t ( ^(x„), t ± e) = Pn (T^(x n ), t) (D.7) 

for any t > 0, which means continuity of the function in (|D.1[) for those t such that T^ n \x n ) £ dT n . 
To deal with the collision configurations, notice that, for all x n € , 

Pn (T t { + l) (x n ), t) = PnpPteJ, t) (D.8) 

for all t > 0, even if this is not true for the initial measure (this property must not be confused with the 
"continuity along trajectories", see Eq. (|4.3[) ). In fact, for n = N Eq. (|E).8|) is again a trivial consequence 
of the Liouville equation, while for n < N, it is easily proved by induction: assuming it for p n +i, 



Pn{T;+{x n ),t) = dx n+ ip n+1 (T} + >(x n ),x n+ i,t) 
•M(T t ( ;>( £ j) 

dx n+1 p n+1 (T t ^ (x n ), x n+ i,t) 

= Pn(T^\x n ),t), (D.9) 

Equation (|D.8|) . together with (|D.7I) . prove the first assertion of the lemma for the correlation functions. 
Coming now to the functions V(T>) and remembering the explicit expression (I3.10p . we observe that 

V(V)(Tt± E ±(x n ), tie) = Jq £ dt\ ■ ■ ■ , where the dots indicate a function that depends only on the states of 
the evolution T> = (T^ e± (x n ) , m, ji , ■ ■ ■ , j m , t x , • • • , t m , p\ , ■ ■ ■ , p m , wi , • • • ,w m ), during the time inter- 
val [0, ti]. Then for any t\ £ (0,t) this function is actually independent on e: we can substitute Tt±\±(x n ) 
in T> with T^{x n ). Thus we obtain V(V) (Tt± e ±(x n ), t ± e) — J* dt\ ■ ■ ■ + f* ±£ dt\ where the first term 
concides with V(V)(T^\x n ), t), and the second term can be bounded, proceeding as after f|3.10[) . with O(e). 
This shows that property (|D.7j) holds also for the function V(T>), while property (|D.8|) is obvious for V(T>), 
so that the claimed continuity property is proved for all x n £ T* and all t > 0. 

Finally, to prove the statement over almost all dT n , it suffices to apply the second part of Lemma [5J □ 
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